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INTRODUCTION 

The promise of filamentary composite materials, whose de- 
velopment may be considered as entering its second generation, 
continues to generate intense interest and applications activ- 
ity. Such interest and activity are well-founded, since they 
are based on the possibility of using relatively brittle mate- 
rials with high modulus, high strength, but low density in 
composites with good durability and high tolerance to damage 
and which, when they do fail, do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the mid 1960* s. Rather limited applications to 
primary structure have been made on operational vehicles , 
however, mainly in a material-substitution mode on military 
aircraft. More extensive experiments are currently underway 
on large passenger airplanes in commercial operation as part 
of NASA's influential ACEE program and in a few military de- 
velopments, such as the AV-8B and X-29 which have not seen 
service use. 

To fulfill the promise of composite materials completely 
requires a strong technology base. NASA and AFOSR recognize 
the present state of the art to be such that to fully exploit 
composites in sophisticated aerospace structures, the tech- 
nology base must be improved. This, in turn, calls for 
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expanding fundamental knowledge and the means by which it 
can be successfully applied in design and manufacture. 

As the technology of composite materials and structures 
moves toward fuller adoption into aerospace structures, some 
of the problems of an earlier era are being solved, others 
which seemed important are being put into perspective as 
relatively minor, and still others unanticipated or put aside 
are emerging as of high priority. The purpose of the RPI 
program as funded by NASA and AFOSR has been to develop crit- 
ical advanced technology in the areas of physical properties, 
structural concepts and analysis, manufacturing, reliability 
and life prediction. 

Our approach to accomplishing these goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing and, where 
the state of the art will be advanced by doing so, extending 
the research effort into simulated service use so that the 
composite structure *s long-term integrity under conditions 
pertinent to such use can be assessed. Inherent in the RPI 
program is the motivation which basic research into the 
structural aspects provides for research at the materials 


level, and vice versa 
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Second, interactions among faculty contributing to pro- 
gram objectives - which is a group wider than that supported 
under the project - is on a day to day basis, regardless of 
organizational lines. Program management is largely at the 
working level, and administrative, scientific and technical 
decisions are made, for the most part, independent of con- 
siderations normally associated with academic departments . 
Involvement of this kind includes - depending on the flow of 
the research - faculty, staff and students from chemistry, 
civil engineering, materials engineering and the department 
of mechanical engineering, aeronautical engineering and 
mechanics . 

Both of these characteristics of the NASA/AFOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

A third aspect is a newly instituted initiative to in- 
crease involvement of NASA's Research Center scientists and 
engineers in the program at RPI and vice versa. This will 
require, first, identification of individual researchers 
within NASA centers whose areas of interest, specialization 
and active investigation are in some way related to those of 
RPI faculty supported under the subject grant. Second, means 
to bring about their interaction will be sought and a program 
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of active interchange encouraged. Important benefits fore- 
seen for this increased communication include a clearer 
window to directions in academia for NASA researchers; oppor- 
tunities to profit from NASA experience, expertise and facil- 
ities for the faculty so involved; and an additional channel 
for cross-fertilization across NASA Research Center missions 
through the campus program. 

Overall program emphasis is on basic, long-term research 
in the following categories: (a) constituent materials, (b) 

composite materials, (c) generic structural elements, (d) 
processing science technology and (e) maintaining long-term 
structural integrity. Emphasis has shifted, and can be ex- 
pected to continue to shift from one time period to another, 
among these areas depending on the states of composite mate- 
rials and structures. Progress in the program will be re- 
ported in the following pages under these headings. Those 
computer methodology developments are also iindertaken which 
both support Rensselaer projects in composite materials and 
structures research in the areas listed above and which also 
represent research with the potential of widely useful re- 
sults in their own right. 

In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scien- 
tists and engineers in a number of pertinent disciplines at 
RPI will interact, both among themselves and with counterpart 
NASA Center researchers, to achieve its goals. Research in 
the basic composition, characteristics and processing science 
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of composite materials and their constituents is balanced 
against the mechanics, conceptual design, fabrication and 
testing of generic structural elements typical of aerospace 
vehicles so as to encourage the discovery of unusual solu- 
tions to present and future problems. In the following 
sections, more detailed descriptions of the progress achieved 
in the various component parts of this comprehensive program 
are presented. 


f 
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II-A CHEMICAL MODIFICATION OF THE SURFACE OF HIGH MODULUS 
CARBON FIBERS 
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I I -A CHEMICAL MODIFICATION OF THE SURFACE OF HIGH MODULUS 
CARBON FIBERS 

Senior Investigator: R. J, Diefendorf 

1 , Introduction 

Over the last two decades, there has been great demand 
for engineering materials with very high specific strength 
and stiffness. Extensive research led to the emergence of 
the families of kevlar, boron, carbon fibers and others. 
Unfortunately, carbon fibers with higher specific modulus 
(>450 GPa) usually have a graphitic surface and are poorly 

rii * 

bonded by most polymeric matrices . The interface between 
high modulus fibers and epoxy matrices accounts for the low 
interfacial shear strength observed in composites comprised 
of these constituents. Consequently, one of the central 
challenges in the development and application of high modulus 
and ultra-high modulus carbon fibers is how to improve their 
composite interlaminar shear strength, while at least main- 
taining composite toughness ^ ; Several approaches to 

modifying either the matrix resin, fiber surface or both have 
met with varying degrees- of success , in that real improve- 
ments may incur high cost, longer processing times or are at 
the expense of composite toughness 

In the past, considerable effort at RPI has been ex- 
pended to study the surface of high modulus fibers , matrix 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 23. 
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resins and epoxy curing agents and to observe how much their 
interfaces differ. In cases where differences do exist, in- 
vestigations were made to try to match these differences and 
to design a more predictable or controlled fiber-matrix in- 
terface. 

The aim of the present sub-study is to describe the 
deposition of selected monomers on the surface of high modu- 
lus carbon fibers. This research, being conducted by gradu- 
ate student Cyprian Uzoh, is part of a continuing effort to 
understand and define the fundamental principles which pro- 
vide the basis for the design of a predictable interfacial 
bond between carbon fibers and polymeric matrices. 

2. Status 

The surface energy of anhydride cured epoxies had been 
characterized and the dominant role of the curing agents and 
meaning of the critical surface energy defined, at the be- 
ginning of this reporting period. 

It was learned in our earlier studies that amines 
preferentially adsorb (or react) on the basal planes of gra- 
phite and dictate the effective value of the surface energy. 
From this work, it became apparent that an effective method 
for the design of a controlled fiber-matrix interface is to 
incorporate or graft to the fiber surface those chemical 
species that could interact with the matrix, with the expec- 
tation that the shear strengths and toughness of the resulting 
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composite could be related to fiber surface chemistry. 

3. Progress During Report Period 

Based on the earlier observations, various chemical 
methods were used to graft functional groups (amine, anhy- 
dride, amide etc.) and monomers on the surface of high modu- 
lus carbon. (Before grafting, the fibers were scoured in a 
50/50 solution of chloroform and chlorobenzene, then with 
chloroform, acetone and, finally, with anhydrous ether, then 
dried at 150® C for 12 hours.) 

The reactivity of the various surface grafts with epoxy 
resins was probed with differential scanning calorimetry 
(DSC) . The type of functional groups on the fiber surface 
was determined with infra-red spectroscopy. The nitrogen 
and hydrogen assay of the various treatments was obtained 
with an elemental analyzer, and, finally, the surface mor- 
phology of the modified fibers was examined with scanning 
electron spectroscopy (SEM) . 

a. Results and Discussion 

Pure graphite is hydrophobic in nature, but this hydro- 
phobicity is decreased if oxygen is present on the surface 
of graphite. Oxidative treatment of carbon fiber surface 
introduces several different polar groups on the fiber sur- 
face. The most abundant species are carboxyl (RCOOH) , phe- 
nolic hydroxy (R-OH) and carbonyl (^C = O) , which are largely 
acidic . The other functional groups known to be 
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present, but not always firmly established, include hydro- 
peroxide, quinones , hydroquinones and lactones. 

Oxygen enrichment of the carbon fiber surface due to oxi- 
dation for brief periods in hydrogen peroxide renders the fi- 
ber surface hydrophilic. Dramatic increases in fiber reac- 
tivity were observed when a peroxide-treated fiber was re- 
acted with diethylenetriamine , in an ultrasonic bath. However, 
the observed weight gains were quite in excess of what would 
have been expected from amine adsorption or reaction with the 
basal and graphitic edge sites on the fiber. The results are 
illustrated in Figure III-A-1 and can be explained by several 
factors, namely: 

a) a small amount of residual peroxide being trap- 
ped between the filaments after it had been 
partially dried, 

b) the wetting of the fiber by the reaction medium 
and 

c) the concentrations of amine in the bath, bath 
temperature and agitation and, to some extent, 
the duration of treatments. 


Solvent-peroxide reaction 

Solute-peroxide reaction 

Solute-fiber surface reaction 

Interaction between by-products 

Interaction between by-products 
and surface 

Figure II-A-1, Schematic of Events Occurring on the Fiber 

During Deposition in a Reactive Medium 
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While the reaction medium, diethyl ether, is known to 
be an inert medium for amines , it is not inert to the resid- 
ual hydrogen peroxide filament’s surface. The 

primary product of the oxidation of aliphatic ether under 

ri31 

mild conditions (30®C-70®C) is hydroperoxide . The pri- 
mary products of hydroperoxide decomposition are esters, al- 

E13 14 151 

cohols, acids and aldehydes ' ' . Apart from the possible 

reactions between the solvent and the peroxide, the nitrogen 
in an amine is in its lowest oxidation state and can be easi- 
ly oxidized by residual peroxide on the fiber surface. Di- 
ethylenetriamine , though a mixed amine, will most likely re- 
act with peroxides to form nitroso compounds, via initially 
formed hydroxyl cimine ^ . A third variation of various feas- 

ible reactions in the vicinity of the oxygen-enriched fiber 

surface is the reaction of amine with some pendant carbonyl 

[171 

group at the edges of graphite planes to form imines 
::rc 0 +RNH2 — = UR + H2O 

The shift of the equilibrium in this reversible reaction can 
be shifted to favor imine formation by the selective removal 
of water. However, the carbonyl group, like other bond 
groups, is hindered on the surface and will only react if an 
amine specie diffuses close enough for them to interact with. 
At lower concentrations of amine, the pickup of amine by the 
fiber is small and not affected by peroxide concentration. 

With increasing amine concentration (>0.25 eg/L) , rapid 
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build-up of oily, gel-like material occurs. Studebaker, in 
a similar experiment using Philblack-0 and Spheron 
concluded that the fixation of nitrogen depends on the oxygen 
containing groups initially present on the surface. The ni- 
trogen was so firmly bonded that only a small fraction could 
be recovered on heating with concentrated sodium hydroxide 
or hydrochloric acid at 300° C, which showed great stability 
of the surface complex formed. 

The presence of ammonium ions in the acidic fiber sur- 
face oxide, for our fiber, is discounted by infrared spec- 
troscopy of surface functional groups. There are sharp re- 
ductions in the peaks of the starting diethylene triamine, 
particularly, N-H stretch and bend (at 2.97y-3.05y and 6.27y) 
and CH stretch and bend (at 3.4y-3.56y and 6.85y), respec- 
tively. Of particular interest is the emergence of a pro- 
nounced peak around 6y. Preliminary wet analysis of fiber 
surface deposit may indicate that the peak is due to the car- 
bonyl group. Even though a comparatively thick deposit is 
not a requirement for the development of a good-coupling be- 
tween modified fibers and matrix, or for the formation of a 
tenacious film, it will be of interest to compare the behav-^ 
ior of composite interface resulting from the reaction of 
amine on a hydrophobic fiber surface to that of an oxygen 
enriched surface. 

The polymerization of diethylene triamine with Nadic 
methyl anhydride (NMA) or with Dodecenylsuccinic anhydride 
(DDSA) to form their corresponding N,N-Disubstituted amide 



17 


on the fiber surface under differenct conditions is shown 
in Figure II-A-2. Fiber weight gain increased with the ef- 
fective concentration [A] [B] of the amine, a base [B] and 
the acid anhydride [A] . However, yields were not as high 
as when anhydrides were absent. As in the case of the de- 
position of amine, rendering the fiber surface hydrophilic, 
by scouring in hydrogen peroxide, increases yield. Perform- 
ing the final deposition (anhydride) in the ultrasonic bath 
reduces fiber weight gain. Finally, the comparatively low 
yield observed for Nadic methyl anhydride on an oxygen- 
enriched fiber surface, is consistent with other results; 
namely, that scouring the fiber surface with hydrogen per- 
oxide renders the surface more acidic. 

Examination of the various treated fibers in a scanning 
electron microscope reveals (as shown in Figures II-A-3a to 
-3e) that, in all instances, lower concentrations of amine 
and anhydride produce fairly uniform deposits. At higher 
concentrations of monomers or weight gains, nonuniform de- 
posits prevail, as represented by the clumps on the fibers 
in Figure II-A-3d. However, when the final deposition is 
performed in an ultrasonic bath, these clumps are dispersed 
(see Figure II-A-3e) . Additional work has been done to re- 
late conditions of treatment to the structure and morphology 
of coating. 



% WEIGHT GAIN 


ON HM CARBON FIBERS 



NOTE: 1) DTA-DDSA; DTA-NMA; dodecenyl succinic 

anhydride and nadic methyl anhydride 
deposition on diethylene amine, res- 
pectively 

2) DTA-DDSA/u; anhydride deposition per 
formed in ultrasonic bath 


3) DTA-NMAH (u,p) ; room temperature depo- 
sition without ultrasonic bath or 
peroxide oxidation 

4) DTA-NMA?^P; no fiber peroxide oxidation 

5) Deposition of only nadic methyl anhy- 
dride on fibers 


Figure II-A-2, Deposition of Amide as a Function of Effec- 
tive Bath Concentration of Amine [B] and 
Anhydride [A] 







Surface of High Modulus Carbon Fibers After 
Scouring with a Solution of Chloroform and 
Chlorobenzene (lOOOx) 



HM Fiber Surface, Preoxidized in Thirty 
Percent Hydrogen Peroxide Before the De 
position of DTA (34.6% gain gel forma 
tion) 
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Figure II-A-3d. 


Fiber Surface After Deposition of DTA 
and DDSA. (fiber weight gain 11.7%) 
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Figure II-A-3e. Fiber Surface After Deposition of DTA and 

DDSA with DDSA Deposition in Ultrasonic 
Bath. (4.8% weight gain) 
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4 , Plans for Upcoming Period 
Considerable effort will be devoted to studying the in- 
fluence of various deposited monomers on the modification of 
fiber-matrix interfacial strength, in high modulus carbon 
fiber-epoxy composites. In particular, the effect of the 
interface on crack propagation is to be studied with single 
fibers embedded in a matrix film. 
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1 1 1- A FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 


1. Introduction 

The deformation and failure behavior of graphite/epoxy 
tubes under biaxial (axial and torsion) loading is being in- 
vestigated. The aim of this research is to increase basic 
understanding of and provide design information for the bi- 
axial response of graphite/epoxy composites to fatigue loads. 

2. Status 

The research activities up to the beginning of this re- 
porting period consisted predominantly of experimental inves- 
tigations of the biaxial fatigue and fracture behavior of 
graphite/epoxy thin-walled tubes. These studies continued 
with tests on [0,±45] graphite/epoxy tubes. We have now 
started literature studies on phenomenological damage accu- 
mulation laws in composites. Based on these studies, a ten- 
tative plan was formulated to study the progression of damage 
by residual strength measurements of prefatigued specimens 
and to establish a multiaxial damage accumulation law. 

3. Progress During Report Period 

II 21 

It is well-known that the Palmgren/Miner ' approach 
to block loadings (one stress or strain amplitude followed by 


Nximbers in brackets in this section refer to the references 
which are listed on page 33. 
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a higher - or lower - amplitude) does not predict a sequence 
effect. In this approach the sequence (high/low or low/high) 
of block loadings has no effect on the predicted life. 

The few results of tests on composites which are avail- 
able indicate that the low/high sequence is more damaging 

rs 41 

than the high/low block loading . 

In establishing phenomenological theories of damage 

[5—81 

accumulation, either residual lifetime ■* or residual 
[4 91 

strength ' are adopted. It was tentatively decided to 
use residual strength as a measure of damage in future 
studies . 

Guided by the Paris law of crack growth and previous 
studies , it is proposed to use 

^=g(D)f(S) (1) 

as an evolution equation for damage. In the above, 

N is the cycle number, 

D is the dimensionless damage, 

S is the dimensionless cimplitude and 

g( ) and f( ) are positive functions which 

are only zero when their argument is 
zero. 

It is postulated that 

— ^f 

D = 0 for a = — — =1 (2) 

^ o 

uts 

and 


D 


for a ^ 0. 


(3) 
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where 

o ^ denotes the fracture strength at any 
time , 

^uts ultimate tensile strength and 

is the dimensionless residual strength. 

For working purposes, it is assumed that the relation- 
ship between and D is given by 

= exp[-h(D)] (4) 

where 

h(D) is a monotonic increasing function of 
argument . 

When special forms for h( ) and g( ) are chosen, the rela- 
tionships proposed in References 3 and 4 can be recovered. 

As an example, the fatigue curve for completely reversed 
loading is expressed as 

S = 1 + r logN (5) 

with r a constant < 0. Further, h(D) = D, and g(D) = 
where a is a value to be chosen and is > than 0. Then the 
residual strength is given by 


( 6 ) 

It can also be shown that Equation (1) predicts block 
sequence independence up to the final stage. In a two-block 
loading sequence the behavior of composites is reproduced by 
(1) together with (2) and (3) . 




1-S 1 1 


0 = exp 


r l-a 

[10 N j 

InS 



32 


A simple extension of Equation (1) to multiaxial loading 
is 

||=g(D)f(S) ■ (7) 

with 

S = F. .S.S. + G.S. (8) 

13 1 3 11 

where 

summation from 1 to 6 is implied over re- 
peated indices, 

and are matrices (vectors) which rep- 
resent the appropriate material anisotropies 
and 

symbol t denotes time. 

For proportional loading with common period the cycle 
dependence proposed in Equation (1) can be recovered. 

Preliminary studies were made using Equations (7) and 
(8) together with results obtained previously . Initial 
correlation is encouraging, and the approach outlined brief- 
ly here will be pursued further. 

4. Plans for Upcoming Period 
The theory described above will continue to be developed 
and checked against data. A test program will be outlined, 
aimed specifically at verifying the above approach and at its 
further development. In addition, biaxial fatigue tests with 
[0,±45J graphite/epoxy specimens under proportional loading 
will be performed. 
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III-B EXPERIMENTAL AND THEORETICAL STUDIES OF MOISTURE AND 
TEMPERATURE EFFECTS ON THE MECHANICAL PROPERTIES OF 
GRAPHITE/EPOXY LAMINATES AND NEAT RESINS 

Senior Investigator : S • S . Sternstein 

1 , Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of the matrix resin. 

Moisture is known to adversely affect the properties of 
both neat epoxy resin and epoxy matrix composites. Inhomo- 
geneous swelling makes a major contribution to the moisture 
degradation of mechanical properties, both in the neat resin 
and the composite- It is postulated that the postcuring 
process can change structure-moisture interactions and par- 
tially alleviate its adverse effects. Both analytical and 
experimental studies have been directed toward understanding 
the physical -mechanical -thermodynamic aspects of this problem. 

2. Status 

Relevant experimental and theoretical findings have been 
summarized in the previous progress report (No. 45, page 35). 
Work reported here has been performed by graduate students 
P. Yang and S. Singh under the direction of Dr. Sternstein. 
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3. Progress During Report Period 

Studies on the effects of moisture and temperature his- 
tories on the mechanical properties of graphite/epoxy compos- 
ite laminates have been concluded during the report period. 

Two theses, one doctoral and one master’s, are in the process 
of being completed. Several journal publications are being 
prepared and will be submitted shortly. A brief summary of 
our findings during the last two years, as generalized in 
the report periods, follows. 

Experimental studies on the effects of postcuring his- 
tory on the sensitivity of moisture uptake and the damage re- 
sulting therefrom have shown that void formation, delamina- 
tion cracks and fiber-matrix interfacial failure are strongly 
dependent on the prior thermal and moisture histories of the 
laminate. A properly postcured laminate displays little or 
no reduction in out-of-plane bending stiffness (hereafter 
referred to as stiffness) when sxibjected to modest moisture 
uptake. Conversely, an as-cured laminate (no postcure) shows 
dramatic reduction in stiffness, which is highly nonlinear 
with moisture uptake, suggesting internal damage. Microscopic 
observations bear out the damage hypothesis and corroborate 
our previous observation that the out-of-phase dynamic me- 
chanical loss modulus (a measure of dissipative work of de- 
formation) increases substantially, when moisture is absorbed 
in an as-cured sample. 

Recognizing the complexity of the curing process and the 
gradients of cure which must exist in any finite size laminate, 
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it is clear that it may be difficult or highly impractical 
to properly postcure real (not laboratory sized) engineering 
structures • . 

The results of detailed observations on the effects of 
moisture uptake, glass transition reduction and reversibility, 
and dynamic loss modulus increase in postcured and as-cured 
laminates will be the subject of the papers to be submitted 
for publication in the near future. 

Our theoretical studies have centered on the application 
of inhomogeneous swelling theory to the moisture induced 
swelling and stress field in epoxy matrix materials. A most 
important result has been obtained in that we have shown that 
a gradient in cross-link density can produce a hydrostatic 
tension field sufficient to cause cavitation, namely, - 
0 ^ = = 16,000 psi . Contrary to what might be expected, 

a gradient in cross-link density is more severe than a spher- 
ical hard particle (a box gradient) and produces much higher 
tensile stresses. Since gel particles and cross-link gradi- 
ents are generally believed to occur in all thermosets, 
these theoretical findings are very promising, in that. they 
form a basis by which environmentally induced (solvent) fail- 
ure due to cavitation can be explained. 

The analysis shows that the stress field resulting from 
an inhomogeneous material structure depends not only on the 
cross-link density at a point but also on the cross-link 
density gradient. Therefore, correlations of solvent uptake. 
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cavitation resistance, dimensional changes (e.g., swelling) 
etc. , with average cross-link density, are invalid. A miss- 
ing internal variable is the cross-link density distribution. 
The chemical and materials characterization of a thermoset, 
which does not evaluate the distribution, can lead to numer- 
ous errors in the evaluation of resin resistance to moisture 
related damage. 

4. Plans for Upcoming Period 
A new apparatus is being designed for the measurement of 
stress and stress state effects on equilibrixim swelling and 
absorption phenomena in multiphase systems including compos- 
ites. Such data are virtually nonexistent in the literature 
and are crucial to the further development of the theory. 
Applications include solvent induced cracking, crazing, fa- 
tigue and delamination. 


5. Current Publications or Presentations by 
Professor Sternstein on this Subject 


"Mechanical Characterizations of Composites" 
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Presented at the Ninth Annual Mechanics of Composites 
Review, Dayton, OH, October 24-26, 1983. 

"Viscoelastic Characterization of Neat Resins and Composites" 

Presented at the Gordon Conference on Composites, 

Santa Barbara, CA, January 16-20, 1984. 

"Mechanical Properties of Composites" 

Presented at the Asilomar Conference on Polymers , 
Asilomar, CA, February 20-22, 1984. 
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III-C NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF COM- 
POSITE FRACTURE 

Senior Investigator: M. S. Shephard 

1, Introduction 

To understand mechanisms of failure in composites it is 
necessary to develop insight into their micromechanical be- 
havior, including interactions between matrix and fibers as 
the load is increased from zero to that corresponding to 
failure* Investigations of these phenomena, either experi- 
mental or numerical, pose serious problems. The purpose of 
this project, being carried out by graduate student Nabil 
Yehia, is to develop a nonlinear finite element analysis 
methodology capable of performing numerical investigations 
into two-dimensional micromechanical failure of composites. 

2. Status 

A generalized program for the two-dimensional analysis 
of static crack growth problems has been developed and is 
being tested. In addition to a detailed examination and de- 
velopment of fracture criteria needed to track cracks in com- 
posites at the micromechanical level, the program employs 
fully automatic finite element mesh generation. This is an 
entirely new approach in program structure that is being 
taken to allow for the automatic tracking of discrete crack 
growth. We found that no existing form . of the maximum strain 
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energy density criterion (i.e. , T-criterion) could be di- 
rectly used for crack propagation tracking, and a modified 
T-criterion was postulated. 

3. Progress During Report Period 

During the last reporting period, our effort concentrated 
on (a) completion of the modified T-criterion, (b) investiga- 
tion of the effect of the singular element size of fracture 
criteria and (c) integration of the component parts of the 
automatic crack propagation system. The system combines fi- 
nite element analysis techniques, fully automatic mesh gener- 
ation and crack propagation criteria into a single program 
capable of analyzing general two-dimensional crack propaga- 
tion. 

The modifications of the geometry data base required to 
account for different crack initiation and propagation situa- 
tions have been developed. Final testing of the system is 
being carried out. 

a. Modified T-Criterion 

The maximum strain energy density criterion, T-criter- 
[ 11 * 

ion , has been found to be a good alternative to the mini- 

[21 

mum strain energy density criterion, S-criterion , for pre- 
dicting crack propagation history. As outlined in the last 
progress report, investigation of the T-criterion has led to 
the modified T-criterion. During this reporting period, 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 61. 
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formulation of the modified T-criterion has been completed, 

including a complete comparison between the S-criterion and 

[31 

the modified T-criterion , For a slant crack problem under 
remote unidirectional load, comparison of the normalized 
fracture loads predicted by the two criterion are shown in 
Figure III-C-1, while the fracture loci are given in Figure 
III-C-2. Note that the modified T-criterion is denoted by R 
in these figures. The S and R criteria yield results which 
are quite close for plane stress problems, while the differ- 
ences are much larger in the case of plane strain. These 
differences arise primarily from the changes in size of plas- 
tic zone associated with plane stress, in one case, and plane 

[31 

strain, in the other . 

b. Singular Element Size Effect 

While the singular element size has been shown to affect 

the numerical evaluation of stress intensity factors, such 

errors appear to be within acceptable engineering accuracy, 

[4 51 

e.g., within 10%, for most cases ' . Errors in stress in- 

tensity factors and in crack propagation angle, obtained in 

[71 

a numerical study of the slant crack problem for practical 
ranges of JJ./a values, are shown in Figures III-C-3 and 4. 

Here too, errors are shown to be within ±10%. On the other 
hand, errors in predicting crack propagation increment for 
this problem go beyond acceptable engineering accuracy for 
small il/a, and the error decreases as Jl/a increases (Figure 
III-C-5) , In the absence of information as to which criterion 
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S-plane stress 
S -plane straiin 



Figure III-C-1, Fracture Load for a Slant Crack Under Remote 

Tension Normalized with Mode I (0 = 90°) 
Fracture Load. (G is the energy release 
rate criterion.) 



S-plane stress 
S-plane strain 

R-plane stress 
R-plane strziin 


UiH 4-1 



for plane strain 
K^/K^ for plane stress 


Figure III-C-2. 


Fracture Loci for a Slant Crack Under Remote 
Tension Normalized with Mode I (3 = 90°) 
Fracture (G is the energy release rate cri- 
terion and a is maximum tangential stress 
criterion. ) 



% Error In stress Intensity factors, Kj and K, 


- 2.0 


Figure III-C-3. Error in Evaluating the Stress Intensity 

Factors in Slant Crack Problems, with (3 = 
Versus the Ratio of Singular Element Side 
Length, I, to Half Crack Length, a. 


/ DCT 

QPDT 


Figure III-C-4. Error in Crack Propagation Angle, 8, Versus 

the Ratio of Va for the Plane Strain, Slant 
Crack Problem, with $ = 40®, Under Remote 
Tension 



% Error in criterion critical value 
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leads to exact results, and since the prediction of crack 
propagation angles and increments are criterion-dependent, 
it is assumed that when the exact stress intensity factors 
are used, the error in predicting crack propagation increment 
and angle is zero. It is noted, therefore, that Figure III- 
C-5 indicates only those errors introduced into the crack 
propagation criteria due to inaccuracies in the stress inten- 
sity factors. A complete study of the £/a effect is given in 
References [6] and [7]- In Figures III-C-3, 4 and 5, QPDT 

r 6 1 

and DCT stand for quarter point displacement techniques^ . 

. . F41 

and displacement correlation techniques ■■ which are two 
techniques used in this study to evaluate the crack tip in- 
tensity factors. 

c. Automatic Crack Tracking System 

The crack propagation system is comprised of two separate 
computer programs. The first is an interactive graphics pro- 
gram used to define the problem to be analyzed. The second 
is the crack propagation analysis program^ which combines 
all aspects of finite element mesh generation and fi- 
nite element analyses. In the problem definition program , 
the geometry of the object and the problem attributes of 
loads, material properties and boiandary conditions are de- 
fined. In addition, fracture- oriented attributes, including 
the existence of crack tips and the values of bond strengths 
at material interfaces, are defined. 

There is no finite element mesh generated during the 
problem definition phase, only the problem geometry in terms 
of material region boundary curves. All the problem attributes 
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are assigned to that geometry. 

The second analysis program automatically carries out 
all facets of the numerical analysis process including; 

1) generating a finite element mesh specifically 

accounting for all crack tips and the applica- 

191 

tion of the problem attributes to that mesh , 

2) carrying out the finite element analysis, 

3) predicting any crack propagation occurring 
during that load increment and 

4) updating the geometry to account for any crack 
growth. 

Since this program accepts a geometric input definition and 
combines mesh generation with the finite element fracture 
analysis procedures, its overall operation and data struc- 
tures are different from those employed in standard finite 
element codes. For example, the combination of the component 
parts of this program is complicated by the fact that differ- 
ent types of algorithms and data interchanges are required 
for each. The modified-quadtree mesh generator makes 
extensive use of recursive algorithms and tree data struc- 
tures, while the finite element and crack propagation portions 
require high precision computation and a large library of 
mathematical functions. To address both these issues, a hy- 
brid programming technique is used. The majority of the mesh 
generation process is written in PASCAL, a language that al- 
lows recursion and the straightforward definition of complex 
data structures, while the majority of the analysis portions 
are written in PORTRAY. 
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The data used by the crack tracking program can be broken 
into three segments, namely; 

1) geometry data: consisting of that information 

that defines the problem geometry and all the 
problem attributes of loads, material proper- 
ties and boundary conditions , 

2) mesh data; consisting of the modif ied-quadtree 
and the finite element mesh tied to it and 

3) analysis data: consisting of all the finite 

element matrices and fracture criterion infor- 
mation. 

Although separate, there are a number of specific ties be- 
tween the three groups of data. The mesh data is created by 
the mesh generator directly from the geometry data, while the 
finite element data is constructed from the mesh data. Fi- 
nally, the geometry data is modified based on the results of 
the finite element analysis and fracture criterion evaluation. 

All program data generation and program operations are 
controlled by the geometry data. During the fracture process 
the model must be repeatedly updated to account for crack 
propagation. These changes are represented in the geometry 
data only, and since all subsequent information is generated 
from that data, all crack initiation and propagation is re- 
flected directly in it, and it is the only data that must be 
updated. A detailed description of the data structure of the 
program can be found in References [8] and 19] . 

Depending on the location of the crack, different modifi- 
cations are made to the geometry data base. Cracks initiating 


50 


Modification 


Initiate crack on 
interior of a 
region 


Before Modification 


After Modification 


Op#rttion 

Number 


Explanation 



1 

Create curve Ci 
( crack flank ) 

2 

Create boundary 61 


Propagate a crack 
on the interior 
of a region 



Extend curve Cl 
( crack flank ) 


Initiate a 
crack on a 
material 
interface 



Divide curve Cl 
into curves 
Cl , C2 ( crack flank ) 
and C3 


Divide curve C2 into 


Propagate 
crack along 
material 
interface 



Propagate 
crack to 
an 

existing 

curve 



4 


curve CZ end t section to 


be merged witti crack flank 


c 

Merge section of 

0 

C2 with Cl 

4 

Divide curve Cl Into 
curves Cl and C3 

3 

Extend curve C2 

6 

Merge boundary B2 
into B1 


Figure III-C-6. Geometry Modification Operations 
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Figure III-C-7. Forty-Five Degree Slant Crack Problem 

(Young's Modulus, fracture toughness, 
thickness and yield stress equal to 1.0, 
Poissons Ratio equal to 0.3.) 
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Figure III-C-8. Starting Mesh for the Slant Crack Problem 




(0.08) (0.15) (0.28) 


Modified T- criterion 
S*crlterlon 


Figure III-C-9. Crack Propagation Path and Corresponding 

Fracture Load (P) Predictions 
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on the interior of a region require the creation of curves 
and boundaries, while cracks initiating at material inter- 
faces require the division of existing curves. Similarly, 
propagation of cracks on the interior of a region require 
extending the curves representing the crack flanks, while 
propagation of cracks along material interfaces requires 
dividing the existing curves and merging a part of the divided 
curves with the ciirve that represents the crack flank. 

The various modifications to the geometry data base can 
be broken down into the following basic operations : 

1) create a curve, 4) divide a curve, 

2) create a boundary, 5) merge two curves and 

3) extend a curve, 6) merge two boundaries. 

Figure III-C-6 illustrates how these operations are combined 
to handle the various situations that arise as cracks initi- 
ate and propagate. Creation of curves and boundaries are 
straightforward additions to the geometry data base. The 
other operations require more complex considerations, but are 
efficiently handled in the hierarchial geometry data structure 
and link lists used to store that information . 

As a preliminary application of the developed system, 
the 45-degree slant crack problem shown in Figure III-C-7 has 
been considered. Figure III-C-8 shows the mesh automatically 
generated for this problem. The crack growth for this prob- 
lem is unstable, as indicated by the negative load increments 

listed in Figure III-C-9 , which are predicted as required for 

. . [3] 

subsequent propagation using either the modified T-criterion 
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f21 

or the S-criterion . The upper crack tip is the first to 
propagate. It hits the exterior boundary, and further prop- 
agation is thus eliminated. The lower left crack tip then 
propagates until complete separation occurs. The mesh of 
the geometry of Figure III-C-9 is shown following crack prop- 
agation in Figure III-C-10 as established by the automatic 
mesh generator; note that the upper crack tip has been elim- 
inated from the mesh. 

As a second example, the plate with two holes, shown 
in Figure III-C-11, was considered. The starting mesh for 
this case is shown in Figure III-C-12. As indicated in this 
figure, one hole is under pressure and the other has a crack 
emanating from it. The crack begins to propagate when the 
pressure, p, equals 2.36. In this case, too, the crack is 
unstable, and the load levels required to cause the next two 
crack propagation increments are only 1.83 and 1.50, res- 
pectively. The crack geometry after four crack propagation 
increments is shown in Figure III-C-13, while the associated 
element mesh is shown in Figure III-C-14. 

4 . Plans for Upcoming Period 

During the next reporting period, effort will concentrate 
on final debugging and testing of the fully automatic meshing 
algorithms integrated with the analysis procedures. In addi- 
tion, attempts will be made to improve the criterion for 
treating cracks on material interfaces, and further micro- 
mechanics analyses will be carried out. 


Figure III-C-10. The Automatically Generated Mesh for the 

Geometry of Figure III-C-9 





Figure III-C-11. Plate with Pressurized Hole (Youngs 

Modulus, fracture toughness/ thickness and 
yield stress equal to 1.0, Poissons Ratio 
equal to 0.3.) 
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Figure III-C-12. Starting Mesh for Plate with Pressurized Hole 
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Figure III-C-13. Geometry for Plate with Pressurized Hole 

After Four Crack Propagation Increments 




Figure III-C-14. Mesh for Plate with Pressuri 2 ed Hole After 

Four Crack Propagation Increments 
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III-D FREE-EDGE FAILURES OF COMPOSITE LAMINATES 

Senior Investigator; T. L. Sham 

1. Introduction 

Analytical study on the free-edge delamination of gra- 
phite/epoxy laminates under static loading is being conducted. 
The objective of this research is to obtain fundamental under- 
standing of the mechanics involved, from the continuum mechan- 
ics point of view, in the process of free-edge delamination of 
composite laminates. Such understanding could serve as a 
basis for control of and design against free-edge delamination, 

2. Status 

The initiation of delamination at free edges has been 
examined using the principles of continuum mechanics. The 
analyses are conducted by idealizing the material response of 
each ply in the laminate as homogeneous , linear elastic and 
anisotropic, A generalized plane strain elastic boundary 
value problem has been formulated to investigate the energy 
variation of such a body. 

3, Progress During Report Period 

Based on the boundary value problem formulated earlier 
and described in the December 1983 progress report, candidate 
global parameters that could serve to characterize the initi- 
ation of free-edge delamination have been investigated. A 
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generalized plane strain finite element program, having the 
capability of analyzing orthotropic layers and oriented at 
arbitrary angles with respect to the loading axis, has been 
developed. The finite element program will be used to assess 
the suitability of the candidate parameters (when used in 
conjunction with a critical-value type of criterion) for 
characterizing the initiation of free-edge de laminations. 

The capability of predicting the effects of stacking sequence, 
thickness effects, etc., on free-edge delamination will be 
the basis for assessing the suitability of these ccuididate 
characterizing parameters. 

4. Plans for Upcoming Period 
Finite element calculations, using the developed pro- 
gram, will be conducted as described in Section 3, above. 



PART IV 

GENERIC STRUCTURAL ELEJIENTS 


IV-A QUANTIFICATION OF SAINT- VENANT ' S PRINCIPLE FOR A 
GENERAL PRISMATIC MEMBER 

IV-B II4PROVED BEAM THEORY FOR ANISOTROPIC MATERIALS 
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IV-A QUANTIFICATION OF SAINT- VENANT ' S PRINCIPLE FOR A GEN- 
ERAL PRISMATIC MEMBER 

Senior Investigator: D. B. Goetschel 

1, Introduction 

Saint-Venant* s principle states that statically equiva- 
lent systems of load on a body will produce nearly identical 
stress fields in regions of the body which are remote from 
the loaded area. Recently, investigators^ ^ have devel- 
oped more precise, mathematical statements of this principle 
to the effect that the difference between the stress fields 
resulting from two statically equivalent load systems shows 
an exponential decay with distance from the loaded area. 

The rate of this decay may be calculated and is of 
critical importance in determining the significance of edge 
effects in structural members. 

Numerous investigators have used finite element 
and finite difference techniques to examine edge and end ef- 
fects for composite structural elements, but these have been 

standard static analyses for particular loadings or boundary 

[71 

conditions. Dong and Goetschel have developed a semiana- 
lytic finite element method for finding the decay rates of 
edge effects in anisotropic layered plates. The purpose of 
the present work is to extend that method to the treatment of 
prismatic members of an arbitrary cross section composed of 
any number of anisotropic elements. This analysis may be 

ic 

Numbers in brackets in this section refer to the references 
which are listed on page 84. 


68 


used to find general solutions to end-effect problems for 
any loading. 

In order to accomplish this, we use the method discussed 
in Reference [7] but perform a two-dimensional analysis 
rather than the one-dimensional analysis used there. We 
adopt the displacement components at the end of the prismatic 
member as the dependent variables. Two-dimensional finite 
element polynomial interpolation functions are used to model 
the behavior through the cross section of the member, and 
exponential behavior is assumed in the direction along the 
axis of the prismatic member. As in Reference [7] , a second- 
order eigenvalue problem results in which the eigenvalues are 
the decay rates and the eigenvectors are the nodal displace- 
ments. Each eigensolution represents one possible mode of 
deformation corresponding to a specific self-equilibrating 
load system which must be determined from the displacement 
field. 

If we use N degrees of freedom (d.o.f.) in the finite 
element scheme, then we have an N x N system of equations 
with N eigensolutions. Having found a number of such solu- 
tions, we may perform a modal superposition analysis to con- 
struct approximate solutions corresponding to any self-equili- 
brating load system. This is analogous to the modal superpo- 
sition analyses used in structural dynamics problems, whereas, 
the methods in References [14-16] are analogous in the same 
way to time-step analyses. 
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2. Status 

Previous to this reporting period, the method described 
had been developed, the governing equations formulated and 
the Finite Element Method computer program written and used 
to find decay rates for several different geometries. 

As we discussed in some detail in the previous report, 
we choose to formulate the governing equations by assuming a 
displacement field of the form 

= P^exp(-yx) (1) 

where 

Ui , u^ and u^ are the x, y and z components of 
the displacement vector at the point (x,y,z) , 

are polynomial interpolation functions in 
the cross section and 

Y is the characteristic decay rate. 

When these displacement functions are substituted into the 

stress equations of equilibrium 

T. . . + f . = 0 (2) 

31/ D 1 

where 

T.. and f. are the components of the stress 

j X 1 

tensor and the body force vector, respectively, 
a repeated index implies summation over 1 , 2 
and 3, and a comma followed by an index denotes 
partial differentiation with respect to the 
corresponding coordinate. 

Finite element discretization leads to a system of equations 


of the form 
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[KHU} = {F} C3) 

where, if we use N degrees of freedom, 

IK] is the N X N stiffness matrix and 

{U} and {F} are N-element column vectors 
representing the displacement d.o.f and 
the corresponding force d.o.f., respectively. 

The stiffness matrix arises from integration of the sur- 
face tractions over the end and sides of each element and of 
the body force over the volume of each element. It is there- 
fore natural to express the stiffness matrix as 

IK] = EK^] + IK^] + IK^] (4) 

where 

e, s and v denote end, side and volume. 

Since we will only be considering problems in which there is 
no body force in the interior of the elements and there are 
no tractions on the sides, we may write Equation (3) as two 
equations : 

[K Hu} = {F} (5) 

e 

(IK^I + [K^]){U) = {0} (6) 

If we expand Equation (6) and gather like powers of y# we ob- 
tain an equation of the form 

(Y^[K2] + YlK^] + [KqJ){U) = {0} (7) 

which is the second-order eigenproblem we must solve . Since 
the matrices [K 2 ] and [K^] are symmetric and [K^] is antisy- 
metric, the eigenvalues, y, and the eigenvectors, {u} , are in 
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general complex. Note that, from Equation (1) , if y is real 
and positive or complex with positive real part then the 
solution decays with distance from the end, whereas, if y 
has negative real part the solution grows exponentially with 
distance from the end. Since we require bounded solutions in 
the semi-infinite prismatic member, solutions with Key < 0 
are not of interest. If y and {u} are complex, then the com- 
plex conjugates y* and.{U}* also provide a solution and both 
must be included in order that the displacement field in the 
member be real. 

Decay rates were calculated for isotropic and homogene- 
ous graphite-epoxy composite plates and the results compared 

[51 

with exact solutions which exist for these cases. Agree- 
ment was found to be excellent. Decay rates have also been 
calculated for isotropic and composite -laminate I-beams, 
box-beams and rectangular-beams, and for a skin-stringer ar- 
rangement. As one would expect, end effects in thin-walled, 
open cross sections decay much more slowly than those in 
solid cross sections; for example, the slowest decay rates 
(smallest Rey*s, corresponding to the most troublesome end 
effects) were up to thirteen times slower for isotropic I- 
beams than for isotropic rectangular beams. It was also 
found that anisotropy can lower the decay rate by a factor of 
three or four for identical cross sections. 
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3, Progress During Report Period 
Since the computer program has been developed and de- 
bugged by the start of this reporting period, recent effort 
has been directed towards extending the postprocessing capa- 
bilities of the program to calculate the displacement and 
stress fields in the prismatic member from the nodal dis- 
placement vectors (eigenvectors) and decay rates (eigen- 
values) and to perform modal superposition analyses to con- 
struct solutions to problems with specific applied forcing 
functions using the calculated eigensolutions , 

The displacement fields corresponding to the decay rates 
with positive real part have been determined for a homogene- 
ous graphite/epoxy plate two units thick. Four elements were 
used, with three nodes per side in the vertical (z) direction 
and two nodes per side in the horizontal (y) direction, as 
shown in Figure IV-A-1. Although the geometry shown appears 
to be that of a beam, in this case the plane strain behavior 
of a plate has been enforced by requiring the y-displacements 
of the nodes to be zero and the x and z displacements for 
nodes with the same y-coordinate to be the same. Thus, for 
eighteen nodes we have 18 d.o.f., and we can determine eight- 
een different modes of deformation. The deformed shapes have 
been plotted and are shown in Figure IV-A-2 for the lowest 
several decay rates. Note that the displacement fields are 
either symmetric or antisymmetric in z and that as Key in- 
creases the deformation becomes increasingly complicated. 



Figure IV-A-l . Geometry for Homogeneous Graphite/Epoxy Plate 


'vl 

CJ 
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Since y is small for the first few modes and the displacement 
is predominantly in the x-direction, it is difficult to see 
the decay of the deformation with x. However, it is clearly 
visible for the higher modes. 

By using the same geometry shown in Figure IV-A-1 but 
allowing all the modes to move freely in any direction, we 
have the behavior of a rectangular beam two units high by 
one-half unit thick. In this case the eighteen nodes give 
us fifty-four different modes of deformation. The corres- 
ponding deformed shapes of the member have been plotted and 
are shown in Figure IV-A-3 for the lowest several decay rates. 
Note that, included in the possible modes of deformation for 
the beam are those for the plate which are shown in Figure 
IV-A-2. However, the geometry of the beam also permits dis- 
placement in the y-direction; note, for example, the bending 
behavior shown in Figure IV-A-3c. 

We have performed modal superposition analyses for spe- 
cific displacement boundary conditions in the following manner. 
Suppose we have an end condition of the displacement field 

'iiC0,y,z) = u^^^^ (y,z) (8) 

(be) 

We must first evaluate u. ^ at the coordinates of each node 

X 

I3 c 

in order to construct a nodal boiindary condition vector {U ^}, 
and then superpose an appropriate number of the eigensolutions 
to form an approximation to {U }. In general we need not 
consider all possible modes of deformation (i.e., all N 
eigensolutions) , since the highest-order modes (largest y*s) 
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do not contribute significantly to the solution for most end 
conditions of practical interest; however, determining what 
number of modes is "appropriate" in a given situation must be 
based to some extent on experience, which is scarce at this 
time. Thus, if we choose n eigensolutions , we can superpose 
them as follows; 

|... |{U<"’)]{a} (9) 

where {U are the nodal displacement vectors of the i-th 
eigensolution, here constituting the columns of a rectangular 
matrix, and {a} is a column vector of undetermined amplitude 
coefficients. {U^^} is a column vector of dimension N (where 
N is the number of d.o.f. at one station, x = 0), is a 

matrix of dimension N x n, and {a} is a colximn vector of di- 
mension n, so that Equation (9) cannot be solved exactly for 
{a} unless N = n; however, in accordance with the discussion 

above, this is rarely necessary. We can obtain a best R.M.S. 

ri7i 

fit by operating on both sides of Equation (9) with 

H (10) 

where superscript T denotes the transpose. Then we have 

{bV°} = [B'^BHa} (11) 

which may readily be solved for {a}. 

Such a modal superposition has been performed for two 
different boundary conditions for the T-section shown in 
Figure IV-A-4a. There are twenty nodes in the isotropic. 
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four-element member, and the bottom two nodes are fixed, so 
that this situation could be thought of as representing the 
symmetric modes of deformation in an I-beam, as well. There 
are 54 degrees-of- freedom. Figures IV-A-4b and 4c show the dis- 
placement boundary conditions imposed in two cases. In both 
cases, the u^-displacement (the x-direction) is zero. In 
Figure IV-A-4b, the flange has been deformed to follow a 
parabolic curve, and the web has been compressed slightly. 

This is a somewhat simpler situation than that shown in Fig- 
ure IV-A-4C, in which the edges of the flange have been 
clamped to form a compound bend, so we would expect to see 
higher-order (larger y) modes play a more important role for 
4c than 4b. 

Some of the modes of deformation we have used in the 
modal superposition analysis are shown in Figures IV-A-5 and 
6. Note that since both deformations are symmetric about 
the plane y = 0, the modes shown in Figure IV-A-6, which are 
antisymmetric about this plane, should not contribute to the 
solutions. In Table IV-A-1 we list the superposition coef- 
ficients [the a's in Eq. (11)] calculated for the modal super- 
position for the simpler boundary condition in Figure IV-A-4b 
using - successively - nineteen modes, thirteen modes, nine 
modes and three modes. Note that the first important mode is 
the fourth. The results of these superpositions plotted 
along with the applied boundary condition are shown in Figure 
IV-A-7. It can be seen that with three modes the agreement 





TABLE IV-A-1 


00 

to 


AMPLITUDE COEFFICIENTS FOR MODAL SUPERPOSITION USING THE END CONDITION IN FIGURE IV-A-5b 


No. 

. . .._T 

19 Modes 

13 Modes 

9 Modes 

1. 

.18421 

-.00593 

-.00104 

-.000397 

o 

.5850 









.9466 




w • 




4. 

1.0215 + .19551 

-1.0989 - 1.05011 

-1.0557- .99811 

-1.008- 1.2111 

5. 

1.0215 - .19551 

-1.0989 + 1.05011 

-1.0557 + .99811 

-1.008+ 1.2111 

6. 

1.0461 + .46941 

.3046 - 1.16251 

.2596 - 1.05911 

.2937-1.1121 

7. 

1.0461 - .46941 

.3046 + 1.16251 

.2596 + 1.05911 

.2937+ 1.1121 

8. 

1.0649+ .59571 

.2613 +.37241 

.2258 +.32221 

.1872+ .32741 

9. 

10. 

1.0649 - .59571 
1.4398 

.2613 -.37241 

.2258 - .32221 

.1872 -.32 741 




11. 

1.8960 - .77131 




12. 

1.8960 + .77131 




13. 

1.9539 

-.003848 

.2711 


14. 

2.9046- .92891 

.0735 - .01031 

— 


15. 

2.9046 + .92891 

.0735 + .01031 



16. 

2.9389- 1.0391 

-.00401 + .07981 



17. 

2.9389 + 1.0391 

-.00401 - .07981 



18. 

2.9640 + 2.2761 




19. 

2.9640 - 2.2761 





3 Modes 
.0224 
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is poor and with nine inodes the agreement is substantially 
improved. With both thirteen and nineteen modes the agree- 
ment is excellent, and the two are virtually indistinguish- 
able. This is because the thirteenth mode is the last im- 
portant one; the amplitudes for the fourteenth through seven- 
teenth modes are several orders of magnitude smaller than 
the thirteenth. Thus, in this case, the "appropriate number" 
of modes we spoke of above is thirteen. 

In Table IV-A-2 are listed the amplitudes for the more 
complicated boundary condition of Figure IV-A-4c. Since the 
boundary conditions are not drastically different from one 
another, this table shows the same trends as Table IV-A-1, 
but it can be seen that the higher-order modes play a some- 
what more significant role in this case. In Figure IV-A-8 
are shown the results of these superpositions along with the 
applied boundary condition. The same remarks made about 
Figure IV-A-7 apply as well to Figure IV-A-8. 

4. Plans for Upcoming Period 
The program will be extended so that traction boundary 
conditions may be specified for modal superposition, as well 
as the displacement boundary conditions considered here. 

5. References 
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TABLE IV-A-2 


AMPLITUDE COEFFICIENTS FOR MODAL SUPERPOSITION USING THE END CONDITION IN FIGURE IV-A-5c 


No. 

Y 

19 Modes 

13 Modes 

9 Modes 

1. 

.18421 

.001875 

-.001214 

.00174 

2. 

.5850 




3. 

.9466 




4. 

1.0215 + .19551 

-1.032 - .74511 

-1.114 - .80581 

-.8919 - 1.7871 

5. 
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7. 
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8. 
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.4865 + 1.8921 
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.4575 - .22821 

.2789 - .25211 

10. 

1.4398 


— 


11. 

1.8960 - .77131 




12. 

1.8960+ .77131 


— 


13. 

1.9539 

1.494 

1.2520 


14. 

2.9046 - .92891 

-.00094 + .04491 


1 

15. 

2.9046+ .92891 

-.00094 - .04491 



16. 

2.9389 - 1.0391 

-.0526 - .05391 



17. 

2.9389+ 1.0391 
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* 

CO 

2.9640 + 2.2761 




19. 

2.9640- 2.2761 





3 Modes 


.0302 
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IV-B IMPROVED BEAM THEORY FOR ANISOTROPIC MATERIALS 

Senior Investigator : 0 . Bauchau 


1, Introduction 

Timoshenko beam theory is widely used by structural 
designers as a first approximation in numerous structural 
applications. For solid cross-sectional beams made out of 
isotropic material, this theory gives accurate predictions 
for aspect ratios L/h >5 (L is the span of the beam and h 
its height) . Because it is based on the assumption that 
cross sections remain plane after deformation, the theory 
predicts a linear distribution of axial strains. However, 
additional axial strains are induced by warping incompati- 
bilities generated either by specific loading or boundary 
conditions, or by the occurrence of nonuniform torsion or 
bending. According to Saint-Venant ' s Principle these addi- 
tional strains tend to decay away from the perturbation that 
created them. The decay length, 6, is defined as the dis- 
tance it takes for these strains to decay to a negligible 
value. For solid cross-sectional beams made out of isotro- 
pic materials, this decay length is of the order of the 
height of the beam, i.e., 6 — h. A physical meaning of the 
validity range of the theory is found by rewriting it as 
L/6 > 5, stating that the span of the beam must be large com- 
pared to the perturbed zone , 6 . 
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In the case of thin-walled box beams, it was recognized 

early that beam theory provides only a poor representation 

because warping and shear lag effects are significant for 

those structures. More recently, with the increasing use of 

composite materials, Saint-Venant' s principle has been inves- 

[1 2 ] * 

tigated theoretically for anisotropic elasticity ' . The 

decay length, 6, was found to be proportional to /E/G (E is 

the longitudinal Young's modulus , and G is the shear modulus). 

This means that for highly anisotropic materials the decay 

length can be much larger than for isotropic materials. 

[3] 

Goetschel calculated these decay lengths for different 
beam configurations made of unidirectional graphite/epoxy. 

For a thin-walled rectangular cross section, 6 was found to 
be 7h, and for an I cross section, 6 = 15h. Beam theory 
can still be applied if L/6 > 5, i.e. L/h > 35 and 75, res- 
pectively. Needless to say, however, such a restriction makes 
beam theory inapplicable to a great many practical structures. 

The purpose of this work is to investigate a more appro- 
priate representation of the behavior of beams made of aniso- 
. tropic material. 


2. Status 

This work was initiated at the beginning of the report- 
ing period, and concentrated on the analysis of thin-walled, 
closed cross-sectional beams made of anisotropic material. 


Numbers in brackets in this section refer to the references 
which are listed on page 99. 
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The analysis is based on the sole assumption that each cross 
section is infinitely rigid in its own plane, but free to 
warp out-of-plane. This assumption has two important impli- 
cations; 1) The in-plane displacements of the cross section 
are fully represented by three rigid body modes, respective- 
ly, i.e. two translations and the rotation of the section. 

2) Any applied transverse load only induces membrane stresses 
in the structure, specifically an axial stress flow, n, and 
a shear stress flow, q. For thin-walled beams, these two 
stress flows are uniform across the thickness and the other 
stress components are assumed negligible. 

Based on these assumptions, a set of orthonormal eigen- 
warping functions can be derived from the principle of minimum 
potential energy. The solution of the problem can then be 
obtained by expanding the axial displacement or axial stress 
distributions as a series of these eigenwarpings . 

3. Progress During Report Period 

During the report period two tasks were completed: 

1) Closed form solutions of the problem have been 
derived using two alternative approaches. In 
the first approach (the "Improved Bernoulli") , 
the axial displacements are expanded as a series 
of eigenwarpings, and the principle of minimum 
potential energy is used to derive the governing 
equations. This solution is an assximed displace- 
ment approach and yields a lower bound of the 
strain energy. In the second approach (the "Im- 
proved St. -Venant"), the axial stresses are 
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expanded as a series of eigenwarpings and the 
Reissener energy principle is used to find the 
solution. Since an equilibrating stress field 
is assumed, the solution yields an upper bound 
of the strain energy, 

2) A finite element solution was also derived us- 
ing the same displacement field as in the "Im- 
proved Bernoulli" approach. The 4-segment 
beam model has the following degrees of freedom 
at each node; 3 displacements, 3 rotations and 
1 degree of freedom per eigenwarping to be con- 
sidered. 

A specific example will be described here to illustrate 
these different approaches. Figure IV-B-1 depicts the thin- 
walled, rectangular cross-sectional beam to be analyzed, and 
Figure IV-B-2 shows the few first eigenwarpings for this 
section. The transverse displacement, v, of the beam under 


a uniform transverse load, p^ , and the twist, 0, under a uni- 


form torque, m , have been calculated using the different 

° V . . 

. tlD 

approaches. Values of the tip deflection, and the tip 


twist, 
remarks : 


6 tip 


mpL 


Po^ 

, listed in Table IV-B-1, call for the following 


-warping effects can be very significant, specially 
in the torsional analysis; 

-all methods have an excellent convergence rate; 

An accurate solution in bending can be obtained 
with 2 eigenwarpings in the Bernoulli approach 
and only 1 eigenwarping in the St.-Venant approach. 
In torsion, convergence requires 3 and 2 eigen- 
warpings, respectively. 
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40 mm 



Figure IV-B-1. The Thin-Walled, Rectangular Cross- 

Sectional Beam Analyzed in the 
Numerical Example 






TABLE IV-B-1 


NUMERICAL RESULTS FOR THE DIFFERENT APPROACHES 


Tip Deflection = [lO^^m] 

Tip Twist = 


P L** 


m L** 

o 


o 

Bernoulli St.~Venant F.E.M, 

Bernoulli St. 

-Venant F.E.M. 


Basic 

Solution 

4.024 

4.860 

4.02 

91.7 

143.3 

91.7 

1st Term 

0.446 

-0.160 

NA 

20.7 

-9.9 

NA 

2nd Term 

0.104 

-0.018 

NA 

10.1 

-2.1 

NA 

3rd Term 

0.028 

-0.003 

NA 

3.6 

-0.5 

NA 

4th Term 

0.002 

-0.000 

NA 

0.7 

-0.1 

NA 

TOTAL 

4.604 

4.667 

4.597 

126.7 

130.8 

126.3 
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-as expected, the Bernoulli solution and the fi- 
nite element approach give similar results, 
since they are based on the same assumed dis- 
placement field. 

Figure IV-B-3 shows the distribution of axial stress, 

— and shear stress, — ^ — , in the upper face of the beam 
. PoL^ 

under a uniform loading, p^ . Instead of the uniform axial 
stress distribution predicted by the basic theory, a large 
shear lag effect is observed. Figure IV-B-4 depicts the 
stress distribution when the loading is a laniformly distri- 
buted torque, ra^. Here the uniform shear stress distribu- 
tion predicted by the basic theory is significantly altered 
and large axial stresses are generated near the corners of 
the section. 


4 . Plans for Upcoming Period 
The main task to be completed in the upcoming period is 
an experimental verification of the improved beam theory pre- 
sented here. Test specimens will be made of two rectangular 
graphite/epoxy facings joined by two aluminum channels to 
fabricate a box beam. Displacements and strains will be 
monitored during a standard static test. A comparison with 
calculated values should provide insights into the accuracy 
and limitations of the different methods. 

Specimens presenting bending-twisting coupling terms 
will be tested, too, to assess the influence of these coup- 
lings on the behavior of the beam. The analytical model will 
be updated to allow these coupling effects to be accounted for. 
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V-A VARIATION OF RESIN PROPERTIES THROUGH THE THICKNESS OF 
CURED SAMPLES 

Senior Investigator: B. Wunderlich 

1 « Introduction 

a, Resin Systems and Reactions 

A commonly used aerospace epoxy matrix is N,N,N',N'- 
tetraglycidyl-4 ,4 ‘-diaminodiphenylmethane (TGDDM) . This is 
also called N,N,N' ,N' -tetraglycidyl-4 ,4 * -methylenebisbenzene- 
amine. It is cured with 4 ,4 * -diaminodiphenylsulphone (DDS) . 
TGDDM, an amorphous, tetraf unctional epoxide monomer, is 
liquid at room temperature, whereas, DDS has a melting tern- 

ni * 

perature of 178°C . The structures of both components 

are given in Table V-A-1. 

Narmco 5208, Fiberite 934 and Hercules 3501-6 are com- 
mercially available TGDDM/DDS resin systems, available as 

prepregs, with the latter two also containing a boron tri- 

\2 31 

f luoride-ethylamine catalyst ' . Commercial TGDDM is sold 

as Araldite MY- 720. The major impurity found in MY-720 is 

the trifunctional (instead of the tetraf unctional) epoxide 
[41 

monomer . The reported formulation for Fiberite 934 is 

64 percent by weight TGDDM (MY-720), 25 percent DDS (Eporal, 

Ciba-Geigy) , 11 percent diglycidyl orthophthalate epoxide 

(DGOP) (GLY-Cel-A-100 , Celanese) and 0.4 percent boron tri- 

(51 

f luoride/ethylamine (Harshaw) . Ciba-Geigy recommends a for- 
mulation of 69.4 percent-by-weight m 720 to 30.6 percent DDS^^^. 

* 

Niombers in brackets in this section refer to the references 
which are listed on page 125. 


TABLE V-A“l 
CHEMICAL FORMULAE 


N,N,N' ,N*-Tetraglycidyl-4 ,4 *-Diaminodiphenyl Methane (TGDDM) 


CHz-CH-CHz 




N 


O 


4 ,4 ' “Diaminodiphenylsulphone (DDS) 
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Curing of TGDDM/DDS systems involves two types of cross- 
linking reactions: epoxide-epoxide and epoxide-amine. Though 

an amine hardener (DDS) is present, epoxide groups will react 

F4 71 . . 

with each other ' . The epoxide-epoxide reaction is prob- 

ably similar to the polymerization of polyethylene oxide. In 
the presence of small amounts of water, epoxides will hydro- 
lyze to glycols , Other active-hydrogen compounds , such as 

[8 91 

amines, will react similarly ' Both acid and base cata- 

lyzed reactions are listed in Tables V-A-2 and 3. The epox- 

[141 

ide-amine reaction is given in Table V-A-4 

b. Mobility in the Epoxy Matrices 

It is believed that the epoxy cross-linked networks form 
into a more or less heterogeneous structure. That is, par- 
ticles of high cross-link density form and grow in the matrix 
of uncross-linked or sparsely cross-linked material. As the 
cross-linking reaction advances, these gel particles become 
interconnected into a global network, but the intrinsic heter- 
ogeneity may remain. Further, the gradual decrease in mobil- 
ity of reacting molecules causes the network built at the end 
of the reaction to be considered lower in cross-link density 
than that which was formed initially. 

The two different cross-link densities should have two 
glass transition temperatures. 

Using differential scanning calorimetry CDSC ) , measure- 
ments of heat capacity can be made to provide a basis for 
study of the glass transition region. Such heat capacity 



TABLE V-A-2 

REACTION EQUATIONS FOR THE POLYMERIZATION WITH ACID 


CRH0CH-CH2 H'" CRHoCH -CH2 

H 


CRH^CHCHjOH 


CRH2CHCH2OH 

+ 

• m 

CRH2CH -CH2 


CRH2CHCH20H 

I 

0 

1 

CRH2CHCH2 


etc. 


Oxide 
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TABLE V-A-3 


REACTION EQUATIONS FOR THE POLYMERIZATION WITH BASE 


Initiation 


— icT^ HOCH 2 CHCH 2 R 
O" 


Propagation 


CH2R CH2R 
I I ^ 

HOCH2C-OCH2C-O 


Olid# 


CH2R CH2R CH2R 

I ( I 

H0CH2C-C0CH2O-0CH2C-0- 


Termination 


CH2R 

H(0CH2C)„-0H + OH- 

H 
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TABLE V-A-4 

CURE REACTION WITH DIAMINE 


Mechanism 


R-NH2+ R-HG - CHo 

\ / 

o 


5^ 

R-NHo 

I ^ 

R-HC - CH. 


R- NH2 

I 

RHC - CH2 

I 

0 “ 


R-NH 

I 

RHC - CH2 
I 

OH 


Reaction 


O 

/ \ 

CH2~CH— CH 2 — + H 2 N 

R 

I 

H2N 


OH 
I I 

CH2--CH-CH2-N 
^ I 
R 
I 


OH 

CH2-CH-CH2'^ 


OH 

I 


N-CH2-CH-CH2^ 


— CH2- CH-CH2 


OH 
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measurements on uncross-linked, partially cross-linked and fully 
cross-linked epoxy resins provide insight into basic mobility 
changes of the molecules. First, the value of the transition 
temperature gives the temperature range where microbronian 
motion starts. Then, knowledge of heat capacity changes at 
the glass transition yields an estimate of the number and 
type mobile units in the molecule. Coupled with information 
about the coefficients of expansion, compressibility and 
group and skeletal vibrations, this can provide a detailed 
understanding of the behavior of the epoxy resin molecule. 

On this basis, we believe that study of changes in heat 
capacity with the advancement of cross-linking provides a 
means for characterizing polymeric networks - 

Although there is a large body of literature concerning 
cross-linked epoxy resins , there appears to be no data avail- 
able on heat capacities of these epoxy materials as such. 

Studies which have been made are limited to reactions which 
occur when these formulations are cross-linked to mechanical, 
thermal and other applications. 

c. Glass Transition and the Hysteresis Phenomenon 

The glass transition, T^, is time dependent. The rate 
at which a liquid is cooled determines the position of T^ and, 
thus, the physical properties of the glass formed. A full 
thermal analysis must consider both the cooling rate Cor the 
thermal history in case of more involved thermal pretreatment) 
and the heating rate through the transition region. The dif- 
ferent cooling and heating rates may cause a hysteresis 
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effect. Fast cooling and then slow heating of the amorphous 
solid through its glass transition region should cause an 
exothermic hysteresis , while slow cooling and then fast heat- 
ing results in an endothermic hysteresis , Both are in- 
dicated in Figures V-A-1 and 2 , respectively. The endother- 
mic hysteresis is usually more easily produced and analyzed 
quantitatively. The hysteresis enthalpy, is the area 

under the hysteresis peak such that a cyclic integral over 
cooling and heating curves (dJCpdT) is zero and satisfies the 
first law of thermodynamics. 

In cross-linked systems the hysteresis peak seems to de- 
crease with increasing cross-link density for the same thermal 
treatment. Decreasing of the hysteresis peak has been found 

in polystyrene with increasing divinyl benzene cross-link 
[121 

density . Similarly, no hysteresis was noted for certain 
[131 

epoxy matrices . A question to be resolved is if, in ef- 
fect, for these samples the glass transition region is tem- 
perature independent within the chosen time scale or if the 
often-observed broadening of the transition region makes ex- 
perimental observation of the hysteresis impossible. 

The width of the glass transition region is usually 
stated in terms of either the difference between the end and 

beginning of the glass transition (T_ and T ), AT_, or be- 

ht 13 z 

tween the extrapolated end and beginning (T 2 and T^) , AT^. 

Figure V-A-3 illustrates the characteristic temperatures of 

the glass transition region, including T the point of half- 

9 
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Figure V-A-1. Exothermic Hysteresis (schematic). 
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Figure V-A-3. Characteristic Temperatures in the Glass Transition Region. 
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vitrification. The end and the beginning of the glass trans- 
ition region are usually least accurate, but indicate the 
often substantial broadening of the glass transition which 
affects the performance of the material. In this study, AT^ 
will be used for the measurement of the broadness of the glass 
transition region with occasional reference to the wider 
range, 


2. Status 

It is the purpose of this work to gain knowledge of the 
glassy materials used as matrices in composites and to study 
the homogeneity resulting from the curing process. We intend 
to link the glass transition quantitatively with the presence 
of a given material. During the last report period, this led 
to the study of the hysteresis phenomenon at various cure 
states. The hysteresis was found to be strongly affected by 
cross-linking. Further testing of these novel approaches is 
planned using equipment and techniques developed in RPI*s 
ATHAS (Advanced Thermal ^alysi^) laboratory. 

As will be described in the next section, the data in 
question will require higher sensitivity calorimetry than is 
currently available. In preparation for this work, we de- 
veloped a better means of measuring heat capacity, leading, 
in some cases, to a precision of ±0.1 percent. A report 
of this work has been prepared for publication in the Journal 
of Thermal Analysis. 
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This project has had the partial support of postdoctoral 
fellows Grebowicz, Suzuki and Menczel, and research assistant 
Mr. Judovits. The project was initiated on May 1, 1983 by 
Mr. Gjaja, and a full report of the work by Gjaja will be 
forthcoming shortly as a doctoral thesis. Dr. Grebowicz and 
Mr. Judovits will continue the work. 

3. Progress During Report Period 

Epoxy resins containing various amounts of hardener 
(TGDDM/DDS system) were cured in a muffle furnace at 473 K 
for seven hours. The glass transition temperature, T^, 
versus weight minus percent of hardener in the epoxy resin 
were measured (see Figure V-A-4) . A limit was rapidly reached 
in T at only two percent hardener. Thus, the glass transi- 
tion of the fully cured epoxy- amine matrix seems not much 
different from the epoxide-epoxide cure. 

To study this catalyst-free curing further, the ver- 
sus cure-time for the epoxide-epoxide reaction was also 
studied. MY 720 was cured by itself in an oil bath at 473 K 
for different lengths of time (Figure V-A-5) . The T^ was 
found to increase exponentially with the cure time, and a 
maximum T of about 450 K was reached after eleven hours - 

g 

The reaction was found to be inhibited by running the sample 
under argon. 

The decrease of the hysteresis with increasing cross- 
link-density was investigated on divinyl benzene cross-linked 
polystyrene to get quantitative information on the structure. 
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Figure V-A 


4 . Glass Transition Temperatures as a Function of Hardener Concentration 
(curing of TGDDM/DDS for 7 hours at 473° K) 
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A relationship between the broadness of the glass transition 
region, AT, and the disappearance of the hysteresis peak, 

AHj^, was noted (see Figure V-A-6) , A similar behavior was 
found for MY 720 cured without catalyst (see Figure V-A-7) . 

In general, with increasing cross-link density, T was found 
to increase, the glass transition region was found to broaden, 
and the hysteresis peak was found to dampen (see Figures V-A-8 
and 9) . At very high cross-link densities, however, the 
trend seems to reverse itself to some degree. That is, with 
increasing cross-linking at high cross-link densities, though 
Tg still increases, the glass transition region begins to 
narrow and the hysteresis peak increases again (see Figure 
V-A-7 and also the squares in Figure V-A-9) . 

With increased cross-link-density an increase in T is 

fi5i 

usually noted . However, a plot of T^ versus percent- by- 
weight of DVB revealed for the first time a minimum of about 
two percent DVB (see Figure V-A-10) • This decrease in T^ 
should be investigated further. It could be due to a plas- 

ri21 

tization effect of the not-fully-reacted monomer, . Pos- 
sibly only one end of the DVB reacted at low DVB concentration 
so that instead of a pure cross-linking effect, which would 

raise T , a branching effect would occur which might decrease 
9 

T . Both these trivial reasons are countered by the extra- 
9 

polation of the low cross-link density curve to zero concen- 
tration which is considerably higher than the normally accept- 
ed T of polystyrene (T = 373 K) . This work connects with 
9 9 
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Figure V-A-6. Change of the Glass Transition Hysteresis Peak Area, 

Breadth of the Glass Transition Region, AT, as a Function of Glass 
Transition, T„, for Cross-Linked Polystyrene. (cooling rate of 
0.5®K/min, hating rate 20® K/min) 
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Figure V-A-7, Change of the Glass Transition Hysteresis Peak Area, 

Breadth of the Glass Transition Region, AT, as a Function of Glass 
Transition, Tg, for the Samples Described in Figure V-A-5 (MY 720) 
(cooling rate^O.5® K/min, heating rate 20® K/min) 
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irregularities found in the analysis of a broad range of PS 
samples analyzed in our laboratory for the lUPAC thermal 
analysis working group. 

4. Plans for Upcoming Period 

The work on changes of with cross-linking will be 
continued. Our initial attempt to link the glass transition 
quantitatively with the presence of a given material for the 
,TGDDM/DDS matrix was complicated by the epoxide-epoxide re- 
action masking the T increase of the fully cured epoxide- 
amine matrix. Since the epoxide-epoxide cured samples are 
more brittle, we will explore how to limit this reaction. 
Initial attempts have proven successful (curing under an in- 
ert atmosphere) . 

After a method has been found to inhibit the epoxide- 
epoxide reaction, the versus hardener characteristics for 
the TGDDM/DDS epoxide-diamine reaction can be studied. 

Finally, the investigations of the decrease of the hys- 
teresis effect and its relationship to the broadness of the 
glass transition will be continued. We will try to simulate 
various microscopic and macroscopic models of hysteresis to 
separate broadening effects from hysteresis. 

Although the original goal of this research has moved 
somewhat into the future, we have gained a better understand- 
ing of cross-linked systems and hope to contribute, after 
solving some of the more fundamental problems, to the applied 
research with correspondingly better preparation. 
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V-B ACCURATE MEASUREMENT OF HEAT CAPACITY BY DIFFERENTIAL 
SCANNING CALORIMETRY 

Senior Investigator: B. Wunderlich 

1 . Introduction 

Experience with high quality heat capacity measurement 
by differential scanning calorimetry is summarized and illus- 
trated, pointing out three major causes of error: 

a) incompatible thermal histories of the sample, 
reference and blank runs; 

b) unstable initial and final isotherms; 

c) incompatible differences between initial and 
final isotherm amplitudes for sample, reference 
and blank runs . 

Considering these problems, it is shown for the case of poly- 
oxymethylene that accuracies in heat capacity of ±0.1 percent 
may be possible. 


2« Status 

Improving the measurement of heat capacity by differ- 
ential scanning calorimetry (DSC) is one of the goals in 
RPI*s laboratory for Advanced Thermal ^alysis (ATHAS) . In 

the past we tested standards for temperature, heat of fusion 

[ 1 ] * 

and heat capacity . Comparisons between various DSC's 

[21 

were made at high and low temperature ranges . Computeri- 
zation was introduced early and compared to standard measure- 

[31 . . 

ments . More recently, this computerization was improved 

* 

Numbers in brackets in this section refer to the references 
which are listed on page 144. 
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. . [41 

to sixteen bit data acquisition , which led over the past 
year to the discovery of some limitations in precision of 
heat capacity measurements. These limitations are described 
here, together with recommendations for avoiding these prob- 
lems. The measurements were made with a Perkin-Elmer DSC-2 , 
but the same effects were seen qualitatively on other instru- 
ments (Perkin-Elmer DSC-4 , Du Pont 990, Mettler 2,000). We 
believe the effect is inherent in differential scanning 
calorimetry. 

In the temperature range above 10® K, adiabatic calori- 
metry was , for a long time , the main technique for precision 
[5] 

heat capacity .. Precision of ±0.1 percent could, with some 
care, be achieved routinely. The main problem in adiabatic 
calorimetry is minimizing and calibrating residual heat 
losses. One good practice is to place the whole calori- 
meter laboratory into a room of constant temperature (±0.5® K) 
to eliminate changes in temperature gradients from within 

the calorimeter to outside it. Differential scanning calori- 
[71 

metry has become precise enough to make acceptable meas- 

[81 

urements of heat capacity since about 1965 *• The much 

larger problem of heat losses in twin calorimetry was by- 
passed by drastically reducing sample mass and increasing 
measuring speed. By now many of the other limiting factors 
of scanning calorimetry, such as data acquisition, sample 
handling and heat flow, have been improved sufficiently that 
a precision of ±0.1 percent by DSC is possible. The remaining 



129 


problem seems to be the handling of transient temperature 
gradients . 

3. Progress During Report Period 

This report deals with problems arising in attempting 

[91 

to measure the heat capacity of liquid poly oxyme thy lene 
over as wide a temperature range as deomposition and crystal- 
lization permits. It includes the work of postdoctoral 
fellow H. Suzuki. Figure V-B-1 provides one illustration. 

Two successive runs, under identical measuring conditions, 
without reloading or repositioning the sample and at the 
same heating rate, gave the results shown. After sufficient 
work to indicate that no change in the sample itself had 
occurred, we had to conclude that instrument effects caused 
curve A to be in error. 

a. Instrumentation 

A Perkin-Elmer DSC-2 was used for all measurements. The 

. [31 

prior Hewlett-Packard calculator used with this instrument 

was replaced with a LMS-IOOO-DSC data station A block 
diagram and data handling flow diagrams of this system have 
been given in Reference [4] . The data station consists of a 
Z-80 based microcomputer with 64 k-bytes of random access 
memory and dual eight-inch floppy disk drives. A multichan- 
nel, 16-bit analog-to-digital converter with differential 
inputs is used for data assimilation. Data collection is 
interrupt-driven, using the DSC programming pulse. There are 




Figure V-B-l. Heat Capacity of Polyoxymethylene as a Function of Temperature 
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at present thirty data points recorded per degree kelvin. 
Typical heat capacity runs make use of fifteen to thirty 
milligrams samples. Aluminum pans of similar mass (±0.1 mg) 
are used for reference and sample. The standard heating 
rate is 10® K/min. All runs are made in a dry nitrogen at- 
mosphere with the calorimeter block being thermostatted at 
about 200® K with a mechanical refrigeration unit. The calor 
imeter is shielded from the room environment, which is thermo 
statted to ±1®K, by the standard dry-box attachment of the 
Perkin-Elmer . 

b. Results 

A large series of DSC- runs (about 100) were made and 
analyzed with respect to their thermal history. 

Study of the runs leading to Figure V-B-1 may serve as 
an example for a large series of similar experiments. It re- 
vealed the three major problems to be discussed: 

a) A sudden change in the (measured) heat flow 
("exotherm") occurs whenever the calorimeter is 

kept at a temperature above the initial isotherm 
for a time before measurement, (The effect 
decreases with increasing temperature differ- 
ence from the isotherm. ) 

b) The heat flow under isothermal condition depends 
on the prior thermal history (for times much 
longer than the obvious approach to steady 
state) . 

c) Serious errors in heat capacity can result from 
coupling sample, reference and blank runs of 
different thermal history (due to different pen 
amplitudes) , 



132 


“ As a first analysis of the "exotherm”, the heat flow 
curves of the sample runs in Figure V-B-1 are shown in Figure 
V-B-2. It is clear from durve A that the sudden deviation 
("exotherm") in heat capacity is due to the deviation in 
heat flow of the sample run in the vicinity of the tempera- 
ture of the prior two-minute isotherm (470® K), The observa- 
tion of such an "exotherm" is not restricted to a particular 
polyoxymethylene sample. The same heat flow pattern can be 
seen even in a reference (sapphire) run and a blank (empty 
pan) run as along as a similar thermal history is given prior 
to the measuring run. Figure V-B-3 depicts examples of the 
"exotherm" observed in sapphire runs. The thermal histories 
given are similar to those for rxin A of Figure V-B-2. A sud- 
den change in heat flow can be seen just above 470° K, the 
temperature at which the calorimeter was kept for two minutes 
prior to the measurement run. Furthermore, this figure il- 
lustrates the time effects on the "exotherm". That is, be- 
fore making a run from 450° K to 500° K, the calorimeter was 
held at the 450° K isotherm for five different periods of 
time from 3 minutes to 10.5 hours. Since the general feature 
of the "exotherm" can be seen in all heat-flow curves, it is 
understood that the prior thermal history does affect the 
actual heat flow, irrespective of the time the calorimeter 
is held at the initial temperature of a run, for times as 
long as ten hours. The reasons why a prior thermal history 
of the calorimeter remains influential for such a long time , 
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Figure V-B-3. Heat Flow Curves Observed on a Sapphire Reference 
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going beyond times needed to decay temperature gradients ^ is 
not clear at present. This is one of the instrumentation 
problems which operators must realize and avoid for quality 
heat capacity measurement. Other similar experiments showed, 
as a matter of course, that the effect decreased with in- 
creasing the temperature difference between the prior iso- 
therm (for the case shown in Figure V-B-1, 470^ K) and the 
initial temperature of a measuring run. 

A problem to be discussed next concerns the different 
isotherm positions at the initial temperature, when that 
temperature is approached from different temperatures. This 
effect is obvious from Figure V-B-2 where different inter- 
cepts with the ordinate are seen for curves A and B. For a 
detailed discussion, computer readings of the isotherm posi- 
tions of all runs used for the experiments of Figure V-B-1 
are listed in columns 2 and 5 of Table V-B-1. Readings at 
455° K for the runs, SAA and SAB correspond to the ordinate 
intercepts in Figure V-B-2, the difference of which amounts 
to 1.1 chart paper divisions. Such a difference in isotherm 
position occurs even when the initial temperature is ap- 
proached in the same direction, but on cooling from differ- 
ent temperatures, i.e., 470° K and 505° K in runs SAA and 
SAB, respectively. More significant differences in the iso- 
therm positions at 455° K can be seen for the runs BL and 
RE, when opposite directions of approach to the initial tem- 
perature are applied in the runs A and B. The different 


136 


TABLE V-B-1 

COMPUTER READINGS OF PEN AMPLITUDES IN RUNS SHOWN IN FIGURE V-B-l 


^mperature 

(Init, ) 



(Final) 

Run ' 

455° K 

470° K 

476° K 

505° K 

BL A 

2960 

2322 

2601 

5539 



(1412) 

(1442) 

[2579] 

B 

1988 

1997 

2382 

5500 



(1045) 

(1081) 

[3512] 

RE A 

2722 

-5897 

-5694 

5258 



(9380) 

(9481) 

[2536] 

B 

1766 

-6230 

-5920 

5171 



(9018) 

(9116) 

[3404] 

SA A 

2062 

-1068 

-1059 

4832 



(3961) 

(4284) 

[2770] 

B 

1319 

-1780 

-1441 

4600 



(4083) 

(4138) 

[3281] 


BL, RE and SA are blank, reference and sample runs, res- 
pectively. 

Values in parentheses are absolute values of pen ampli- 
tudes (amplitude minus baseline, interpolated between 
initial and final isotherms) . 

Values in square brackets are differences between initial 
and final isotherms. 

For run conditions, see legend to Figure V-B-1. 

One chart division (2 mm) - 655 in computer units when 
the DSC settings are as follows; gain = 100, range - 
0,5 m cal/s and recorder range = 200 mV, 
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isotherm positions observed after the attainment of steady 
state did not change with time over periods of more than ten 
hours. In contrast, it is seen from column 5 of Table V^B-1 
that the differences in isotherm positions at the final tem- 
peratures of corresponding runs are negligibly small. This 
means that the slope of the line connecting initial and final 
isotherm positions changes with different approaches to the 
initial isotherm temperature. Hence, this affects the actual 
pen amplitude and, accordingly, the value of the heat capa- 
city. In order to examine the possible error in heat capa- 
city, computer readings of the pen amplitude at two interme- 
diate temperatures, 470° K and 476° K, are listed in Table 
V-B-1 in columns 3 and 4, respectively. 

The absolute values of the pen amplitudes (in parenthe- 
ses) in Table V-B-1 should be invariant for runs A and B. 

Yet, they are not. For example, in the runs of RE A and REB, 
the absolute values of the pen amplitude at 470° K are 9380 
and 9018, respectively. This difference of four percent 
is caused by the different isotherm positions at 455° K in 
runs A and B. Assuming identical weights of the aluminum 
pans used, the heat capacity of the sample is proportional 
to (SA-BL/ (RE-BL) , where SA, BL and RE stand for the respec- 
tive absolute pen amplitudes of sample, blank and reference. 
At 470° K this computation yields 0.3199 for runs A and 
0.3810 for runs B. Accordingly, the heat capacity obtained 
from runs A is seen to be sixteen percent less than that 
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from runs B. Similarly, at 476° K, runs A yield a heat ca- 
pacity value seven percent smaller than do runs B, These 
differences in heat capacity are those which have already 
been noted in Figure V-B-1. Their magnitude is aggravated 
by the relatively small sample mass used for the present il- 
lustration. The volume of the sample pan would permit a 
sample four times larger, and that would cut the errors to 
one-quarter the values shown, but still much too large for 
precision heat capacity. 

In addition to the three major problems just discussed, 
there exist two additional minor effects. One of them is 
the roughness of the recording, as can be seen in Figure V- 
B-1. Those data were recorded when condensation was present 
on the surface of the DSC head. Since the curves become 
much smoother after defrosting and cleaning, it can be said 
that condensation introduces noise into the recording. The 
other effect is illustrated in Figure V-B-4. On slow runs, 
a discontinuous adjustment of the average temperature appears 
in the heat-flow curve every 10° K . At heating rates of 10° 
K/min or faster, the steps in the heat- flow curve are practi- 
cally lost. The steps remain the same at all range settings. 
¥ 

c. Conclusions 

From the above findings and considerations, one may set 
up a guide to better DSC experiments: prior to runs, fixed, 

identical thermal histories must be decided upon for sample, 
reference and blank runs. Preferably, we suggest heating to 




Figure V-B-4, Plot of Heat Plow Versus Temperature 
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the final temperature of a run, to erase any prior history 
and then cooling to the initial isotherm at a given rate. 

Even if this precaution is taken, all results may not be us- 
able due to other inherent, well-known problems encountered 
in attempting to achieve reproducibility among runs (pan and 
lid positioning, gas flow, condensation, sample shifting, 
radiation changes by surface contamination, etc.). Only 
without reloading and repositioning a pan does one achieve 
repeated runs of almost perfect reproducibility. Strictly 
speaking, the same experimental conditions cannot be kept for 
blank, reference and sample runs. What must be managed is to 
minimize the changes in experimental conditions and to as- 
sess the quality of the obtained data from the following 
three point procedure. The first point is a check of whether 
three different runs have the same thermal history. If they 
have, they may be called compatible. The second point re- 
quires the check of initial and final isotherm stability. 

The third point involves a comparison of the actually achieved 
differences between the initial and final isotherm amplitudes, 
at least for reference and sample runs, but preferably for 
all three runs. Only if all three checks are positive can 
runs be considered of good quality. Generally, the baseline 
always tilts to some extent. This slope should remain con- 
stant during a given set of experiments. 

As a final example, data on molten polyoxymethylene are 
presented. This polymer can be heated up to 543° K at rates 
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of 40® K/min only once, since on prolonged exposure to high 
temperatures it is subject to significant deoraposition. Fig- 
ure V-B-5 shows three sets of heat-capacities calculated for 
one sample run, H14D1, coupled with different blank and re- 
ference runs. The isotherm amplitudes of all runs are listed 
in Table V-B-2, and the run conditions are described in the 
footnote of the table. Curve A in Figure V-B-5 shows the re- 
sults obtained from the incompatible set of runs H14D1, RES 
and BL3 (the latter two runs have a different thermal his- 
tory) • The result - namely, that the heat capacity A of the 
melt polyoxymethylene is almost independent of temperature - 
must be in error. Curves B and C, practically indistinguish- 
able in the figure, are obtained by coupling the sample run 
H14D1 with RES and BLl and with REl and BLl, respectively. 

All these runs are compatible with H14D1 and, needless to 
say, have stable isotherms. The run BL5 was not adopted, 
because its difference in isotherm amplitudes is too large, 
for unknown reasons. It is seen from column 5 that the dif- 
ference in isotherm position of RES is closer to that of 
H14D1 than REl. It can be seen from Figure V-B-5, however, 
that the observed ten percent difference in slope of the 
reference run is not sufficient to yield a noticeable differ- 
ence in heat capacity. 

The data obtained under such maximized conditions, which 
may require multiple baseline and reference runs , may push 
the precision of heat capacity runs to ±0.1 percent. This 



TABLE V-B-2 


COMPUTER READINGS OF THE ISOTHERM AMPLITUDES 
OF THE RUNS SHOWN IN FIGURE V-B-5 


(Init • ) 
457° K 

191.4 
30.2 
139.8 

RE 1 143.0 

3 12.4 

5 127.6 

H14D1 100.0 


(Final) Difference 

543° K ° K 


546.2 

354.8 

527.6 

497.4 

666.4 

526.6 

472.2 

329.2 

464.4 

452.0 

425.2 

297.6 

401.0 

301.0 



In all runs except BL 3 and RE 3, the pans were 
loaded into the holder at 400° K at 10° K/min. In 
runs of BL 3 and RE 3, the pan was cooled from 
543° K to 457° K at 10° K/min prior to each run. 
Run conditions are the following: Heating rate = 

40° K/min, range = 0.5 mcal/sec, range on the re- 
corder = 500 mV, gain = 10. On these settings, 
one chart division ^70 in computer units. 

H14D1 is the name of a sample run of polyoxy- 
methylene. 
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Figure V-B-5. Heat Capacity of Polyoxymethylene Melt as a Function of Temperature 
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extreme care and large amount of calibration suggests a need 
for advances in instrumentation which, in the past ten years, 
may have been neglected in favor of the development of com- 
puterization. 


4 . Plans for Upcoming Period 
This report concludes the research activities associ- 
ated with heat capacity measurement using DSC. The tech- 
niques developed here will now be used to ascertain the un- 
forraity of cure through the thickness of epoxy matrix com- 
posite laminates. 
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V-C INITIAL SAILPLANE PROJECT; THE RP~1 

Senior Investigators : F. P. Bundy 

R.J. Diefendorf 
H. Hagerup 

On Septeinber 15, 19 83 the annual static loading tests 
were made on the assembled wing/fuselage structure. Simple 
wing bending and wing bending/torsion tests up to about four 
G's were made, using sand bag loading on the inverted air- 
craft. Bending and torsion deflections were measured and 
compared to those obtained in previous years* tests. It ap- 
pears that the structure is retaining its original strength 
and stiffness very well. 

This aircraft has been disassembled and stored on the 
loft of the Composite Materials Shop Area in the Jonsson 
Engineering Center under ambient, indoor conditions of tem- 
perature and humidity. It will be given its annual strength 
and stiffness tests again about August or September, 1984. 
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V-D SECOND SAILPLANE PROJECT; THE RP’-2 

Senior Investigators: F. P. Bundy 

R. J. Diefendorf 
H. Hagerup 

1, Status 

In the previous report the results and analysis of the 
static strength proof tests of the sing/fuselage assembly 
were presented. The configuration tested was that of the 
glider with rebuilt wing, after failure in an earlier test, 
and redesigned carry-through and fuselage load-transfer 
structure. Except for excessive forward wing-bending deflec- 
tions, which caused about a 3/16th-inch pull-out of the aft 
wing pins from their fuselage receptacles when loaded to 
about six G*s, the modified and rebuilt airframe showed very 
good strength and stiffness characteristics. The redesigned 
and reconstructed wing beam carry-through, which transfers 
bending moments via two major shear pins oriented with their 
axes fore and aft through the webs of the port and starboard 
wing spar-roots, withstood the 6-G tests very well. 

2. Progress During Report Period 
During the present reporting period progress has been 
made on the design fabrication and installation of a tension 
carry-through linkage between the aft parts of the port and 
starboard wings to resist forward bending deflection. In ad- 
dition, towhook/release mechanisms, landing gear and wheel 
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brake asseinbly, compound curved sandwich composite fuselage 
shells and an enclosed trailer to house and transport the 
completed disassembled aircraft, have all undergone various 
stages of design, fabrication, installation and testing* 

Using the female molds constructed during previous re- 
port periods, two sets of left and right fuselage half-shells 
were, fabricated. The first set was composite sandwich with 
the layers being, from outside to inside; medium Kevlar 
cloth, graphite cloth, medium Kevlar cloth, Divinycell foam, 
glass fiber cloth. This layering resulted in a shell with 
reasonably good overall contours but with several undesirable 
irregularities. A second set was made which utilized differ- 
ent layering and better vacuum bagging techniques. The lay- 
ering, outside to inside, was as follows: ±45° glass cloth, 

±45° Kevlar cloth, 0-90° Kevlar cloth, 1/8" Divinycell foam, 
±45° Kevlar cloth. This produced a shell with very few minor 
surface imperfections, very good strength and stiffness, and 
a very satisfactory total weight of only eleven pounds. 

An aft wing tension carry-through linkage was designed 
which was intended to (i) be consistent with the existing 
fuselage frame and wing structures; (ii) provide adequate 
tension strength, (iii) provide easy access for pinning and 
unpinning during assembly or disassembly of the sailplane 
and (iv) be consistent with clearances available during stor- 
age in the trailer. A full-size test model of the portions 
of the aft/inboard part of the wing which are involved was 
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made with the appropriate linkage fastenings to the wing 
skins and ribs. This model was load tested in the Instron 
machine and found to be fully adequate. This system is now 
being installed in the wings and fuselage. The stay-pin ac- 
cess is gained by opening the split flap of the wing to the 
full sixty-degree position. 

The towhook/release mechanisms have been fabricated, 
installed and tested for the magnitudes and directions of 
tow cable tensions specified by the FAA, These tests veri- 
fied the capability of the fuselage structure to transmit 
towline force to the wing roots, as in actual flight; and of 
the release mechanism to release the towline when under full 
load in the various specified directions - 

The fuselage/landing gear system was statically tested 
for the FAA-specified capabilities of four-G up, two-G aft 
and 2/3-G sidewise. The system passed these tests without 
indication of any problems. 

The fuselage/tail boom/erapennage structure was given 
the specified FAA strength tests by rigidly mounting the struc- 
ture by the four wing pin sockets and applying sand bag load- 
ing to the empennage airfoils up to the twelve pounds per 
square foot (psf) appropriate to the top speed allowed by the 
flight envelope for this aircraft. When the vertical stabi- 
lizer/rudder was so loaded the assembly was subjected to com- 
bined bending and torsion. This test was also survived suc- 
cessfully. For the horizontal stabilizer/elevator test, the 
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sandbag loading was symmetrical, left and right, so that only 
bending was involved. The structure passed this test, too, 
without any evidence of weakness. 

An invitation was made by the Soaring Society of America 
(SSA) to RPI to exhibit the RP-2 glider at its annual conven- 
tion, March 1-4, 1984, at Hartford, Connecticut. We accepted, 
and, thus, it became necessary to complete an enclosed trailer 
to transport the sailplane on the highways. A trailer was 
constructed during January and February of 1984 as a mono- 
coque rectangular box with appropriate running gear, hitch, 
lights and internal fixtures to hold the glider component 
parts securely during transport. This trailer is now avail- 
able for convenient and secure storage and transport of the 
sailplane. 

At the SSA Convention, the RP-2 was among more than 
twenty other aircraft of sailplane and powered- sailplane 
categories exhibited. The RP-2 attracted considerable inter- 
est on the part of many dozens of attendees. At the technical 
lecture sessions, a paper was presented (by invitation) on 
the RPI Composite Materials Glider Projects by Dr. F. P. 

Bundy. It was received with lively interest by people such 
as Paul McCready (the Gossamer Condor man-powered airplane) , 
Dick Johnson (eight times U. S. National Sailplane Champion 
and well-known aeronautical engineer) and Paul Schweizer (re- 
tired former President of Schweizer Aircraft Corporation) . 
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3. Plans for Upcoming Period 

Work will now be directed toward completion of the sail- 
plane to flight-readiness condition. The various tasks re- 
maining are: 

a) Completion of the cockpit controls (rudder pedal 
system, control stick trim mechanism, flap notch 
selector) , 

b) application of the fuselage shells, 

c) fabrication of the cockpit rim and canopy frame 
to fit it, 

d) complete installation of the aft wing tension 
linkage , 

e) installation and fitting out of the instrument 
panel and 

f) filling, smoothing and painting. 

Before any flight testing is started some additional 
tests will have to be performed. One will be the application 
of a full six-G simple bending load to the wing-fuselage 
structure, which will incorporate the tension carry- through 
near the aft wing pins. Another test will involve a combined 
wing bend- twist loading, simulating conditions at the bottom- 
right corner of the flight envelope (96 mph, -3.8 G*s, -7000 
Ib/in torque) . The weight and center of gravity location for 
the completed aircraft must then be measured prior to flight 
testing. 

Finally, when the aircraft is completed, static-tested, 
weighed and balanced, flight tests will be conducted at 
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Saratoga County Airport by automobile tow, winch tow and/or 
airplane tow. 


4. Current Publications or Presentations by 
Professor Bundy on this Subject 

"Composite Material Glider Program at RPI" 

Presented at Soaring Society of America Convention, 
Hartford, CT, March 1-4, 1984. 


PART VI 

TECHNICAL INTERCHANGE 
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Interchange of technical information is provided through 
technical meetings, both on- and off-campus. A central 
listing of upcoming meetings is compiled, maintained and 
distributed periodically, in order to assure that a Rensse- 
laer faculty/staff member can participate. The calendar 
for this reporting period is shown in Table VI-1, Table 
VI -2 shows the meetings attended by composites program 
faculty/staf f/students from RPI during the reporting period. 
Some on-campus meetings, with special speakers particularly 
relevant to composites, are listed in Table VI-3, A list of 
composites-related visits by RPI faculty/staf f/students , to 
relevant organizations along with the purpose of each visit, 
is presented in Table VI-4, 

The diversity of the research conducted within this pro- 
gram has increased over the last several years. To insure 
information transfer, once-a-week luncheon programs have 
been held among the faculty and graduate students involved 
(listed in Part VII - Personnel - of this report) . These 
meetings are held during the academic year and are known as 
"Brown Bag Lunches" (BBL*s), since attendees bring their own. 
Each BBL allows an opportunity for graduate students and 
faculty to briefly present plans for, problems encountered 
in and recent results from their individual projects , These 
seminars also are occasions for brief reports on the content 
of off-campus meetings attended by any of the faculty/staff 
participants (as listed in Tables VI-2 and VI-4 of this re- 
port) and for brief administrative reports, usually on the 
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TABLE VI -1 

CALENDAR OF COMPOSITES -RELATED MEETINGS 
(September 30, 1983 through April 30, 1984 


1983 

10/17-19 

10/17-19 

10/24-26 

10/24-26 

11/6-10 

11/7-10 


Aircraft Design Systems and Operations Meeting, Fort 
Worth, TX. Sponsored by AIAA/AHS. 

Society of Rheology Meeting, Knoxville, TN. 

Computers in Aerospace IV, Hartford, CT. Sponsored 
by AIAA. 

9th Annual Mechanics of Composites Review, Dayton, 
OH. Sponsored by AFOSR. 

Recent Advances in Experimental Characterization of 
Composites, Salt Lake City, UT. . Sponsored by SESA. 

4th International Conference on Metal Matrix Compos- 
ite Materials, Bratislava, Czechloslavakia. 


11/13-18 Symposium on Impact of New Computing Systems on 

Computational Mechanics, Boston, MA. Sponsored by 
ASME. 

12/5-9 International Conference on Carbon Fibers, San Jose 
dos Compos , Brazil . 


1984 

1/9-11 

1/9-12 

1/9-13 

1/16-20 

2 / 20-22 

2/27-29 

3/14 

3/8-13 

3/26-30 

4/3-5 

4/3-5 

4/6 


Experimental Mechanics of Fiber Reinforced Composite 
Materials, San Diego, CA. Sponsored by SESA. 

22nd Aerospace Sciences Meeting, Reno, NV. Sponsored 
by AIAA. 

Gordon Research Conference on Polymers, Santa Barbara, 
CA. 

Gordon Research Conference on Composites, Santa Bar- 
bara, CA. 

Asilomar Conference on Polymers, Asilomar, CA. 

Vibration Damping Workshop, Applications to Aerospace 
Structures, Long Beach, CA. Sponsored by Flight Dy- 
namics Lab. , AF Wright Aero. Lab. 

Soaring Society of America, Hartford, CT. 

American Chemical Society, Detroit, MI. 

American Physical Society, St. Louis, MO. 

International Conference on Limits in Fiber Reinforced 
Composites, Liverpool, UK. Sponsored by RAES-PRI. 

29th National Symposium and Exhibition on Materials 
and Processes, Reno, NV. Sponsored by SAMPE. 

3rd International Conference on Computing in Civil 
Engineering, San Diego, CA. 


\ 
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TABLE VI -2 

COMPOSITES-RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(September 30, 1983 through April 30, 1984) 

1983 

10/17-19 Society of Rheology Meeting (Prof, Sternstein) , 
Knoxville, TN. 

Professor Sternstein presented the paper: 
"Mechanical Characterization of Composites". 

10/24-26 9th Annual Mechanics of Composites Review (Prof. 
Bauchau and Prof. Sternstein) , Dayton, OH. 

Professor Sternstein presented the paper: 
"Overview of Composites Research". 


11/7-10 4th Inti. Conference on Metal Matrix Composite 

Materials (Prof, Diefendorf) , Bratislava, Czech- 
oslovakia, 

Professor Diefendorf presented the papers: 

"Consumer Applications of Carbon Fiber Composites", 
"Transverse Composite and Fiber Properties". 

12/5-9 Inti, Conference on Carbon Fibers (Prof. Diefen- 
dorf), San Jose dos Compos, Brazil. 

Professor Diefendorf presented the papers: 

"A Comparison of High Modulus Fibers", "Non-Aero- 
space Applications of Composite Materials". 


1984 

1/16-20 Gordon Conference on Composites (Prof. Sternstein) , 
Santa Barbara, CA, 

Professor Sternstein presented the paper; 
"Viscoelastic Characterization of Neat Resins and 
Composites" . 

2/20—22 Asilomar Conference on Polymers (Prof. Sternstein) , 
Asilomar, CA, 

Professor Sternstein presented the paper: 
"Mechanical Properties of Composites", 

3/1-4 Soaring Society of America (Prof. Bundy) , Hartford, 

CT. 

Professor Bundy presented the paper: 

"Composite Materials Glider Program at RPI, 
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TABLE VI- 


1984 

3/8-13 


3/26-30 


4/6 


2 continued 


American Chemical Society Meeting (Prof. Diefen- 
dorf) (Prof. Wunderlich with Drs. Suzuki and Grebo- 
wicz ) , Detroit, MI. 

Professor Diefendorf presented the paper: 

"Carbon Fibers". 

American Physical Society (Prof. Wunderlich with 
Drs. Suzuki and Grebowicz) , St. Louis, MO. 

Presented reports on RPI*s general heat capacity 
work. 

3rd Inti. Conference on Computing in Civil Engi- 
neering (Prof. Shephard) , San Diego, CA. 

Professor Shephard presented the paper: 

"Finite Element Modelling Developments at RPI". 
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TABLE VI -3 

COMPOSITES^KELATED MEETINGS/TALKS HELD AT RPI 
(September 30, 1983 through April 30, 1984) 


Topic 


"Fundamental Studies of 
Acoustic Emission - The Role 
of Theory in the Design of 
Experiments" 


"Boundary Element Method 
Applications to Small and 
Large Deformation Visco- 
plasticity" 

Discussion of proposed 
research into manufacturing 
methods for advanced compos- 
ites 

"Thermomechanical Analysis 
of the Compression-Molding 
Process for Chopped Fiber 
Composites" 

"Plasticity and Cracking in 
Composite Laminates" 


"Asymptotics and Mechanics" 


"Thermoelasticity of Twin- 
ning" 


"Mechanics of Crystalline 
Plasticity" 


Date Speaker (s) 

1983 

10/25 Dr. Nelson Hsu 

National Bureau of Stan- 
dards, Washington, D.C. 

1984 

2/14 Professor S. Mukherjee 

Dept, of Theoretical and 
Applied Mechanics , Cornell 
University, Ithaca, NY 

2/16 Dr. H. Burte and Dr. V. 
Russo, USAF Materials 
Laboratory, WPAFB, Day- 
ton , OH 

2/28 Dr. D. Caulk 

Engineering Mechanics 
Department, General 
Motors Research Labs . 

3/2 Dr. G. J, Dvorak 

Dept, of Civil Engineer- 
ing, University of Utah, 
Provo , UT 

3/22 Dr. J. D. Cole 

Dept, of Mathematics, 
Rensselaer Polytechnic 
Institute, Troy, NY 

3/23 Prof. J. L. Ericksen 

Depts . of Mathematics and 
Mechanics, University of 
Minnesota, Minneapolis, MN 

Dr. A. Needleman 
Prof, of Engineering, 

Brown University, 
Providence , RI 


4/4 
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TABLE VI -4 

COMPOSITES-RELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculty /Staff/Students 
(Septeinber 30, 1983 through April 30, 1984) 


Visited 

Date 

By 

Purpose 


1983 



Fort Lewis and 
AMMRC, Fort 
Lewis, WA 

10/20 

Prof. R. J. 
Diefendorf 

Presented a talk : 
"Composite Structures 
for Army Applications" 

Saratoga Springs , 
NY 

11/5 

Prof. R. J. 
Diefendorf 

Conducted GE Modern 
Engineering Course; 
"Composite Materials" 

Lockheed-Georgia 
Co., Marietta, 

GA 

11/15 

Prof. R. G. 
Loewy 

A review of company 
progress in composite 
technologies 

NASA Langley 
Research Center, 
Langley, VA 

12/6 

Prof. R. G. 
Loewy 

Attended USAF SAB Meet- 
ing on "The Health of 
Aeronautics " 


1984 



NASA and DOD, 
Austin, TX 

1/15- 

19 

Prof. R. J. 
Diefendorf 

Presented a Summary of 
Materials Panel : "Mate 

rials 2000" 

NASA and DOD, 
Washington , DC 

3/1 

Prof. R. J. 
Diefendorf 

Same as above 

USAF Office of 
Scientific Re- 
search, Bolling 
AFB, VA 

3/15 

Prof. R. G. 
Loewy 

Discussion of RPI re- 
search in composites 
structural analysis 

GE Research and 
Development Cen- 
ter, Schenectady, 
NY 

4/17 

Prof. R. J. 
Diefendorf 

Presented a talk: 

"RP Sailplane Project" 

Naval Surface 
VJeapons Center, 

4/24- 

25 

Prof. R. J. 
Diefendorf 

Presented an Aluminum/ 
Carbon Composite Work- 


White Oaks, CA shop: "Carbon Fiber 

Structures and Proper 
ties" 
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part of one of the Co-Principal Investigators* Off-campus 
visitors, at RPI during a BBL day, are often invited to "sit 
in". Table VI-5 lists a calendar of internal, oral progress 
reports as they were given at BBL*s during this reporting 
period. 

An integral part of the program management is an on- 
site review of research plans and progress, conducted under 
the aegis of the technical monitors at NASA and AFOSR. Such 
a program review was conducted at RPI on December 19 and 20, 
1983. The. agenda for this review is shown in Table VI-6, 
and the attendees from AFOSR; NASA headquarters, Langley Re- 
search Center and Lewis Research Center; and industry (mem- 
bers of the program's Industrial Technical Advisory Panel) 
are shown in Table VI-7. 
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TABLE VI -5 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
BROWN BAG LUNCH (BBL) SCHEDULE 
(September 30, 1983 through April 30, 1984) 


Date Topic 

1983 

9/30 Administrative Report 

Fabrication Technology Experiments 


Static Test Review 


Responsible 

Faculty 


R. Loewy 

F . Bundy 
H. Hager up 
V. Paedelt 

H. Scarton 


10/7 Administrative Report 

Resin-Fiber Interface Research 
Composites Fatigue Research 

10/14 Administrative Report 

Matrix Characterization and Environ- 
mental Effects 

Generic Structural Components 


J. Diefendorf 
J. Diefendorf 
E , Krempl 

J. Diefendorf 

S. Sternstein 
D. Goetschel 


10/21 Administrative Report 


J. Diefendorf 


Computer Aided Design and Analysis 


M. Shephard 


Fabrication Technology Experiments 


F . Bundy 
H. Hager up 
V. Paedelt 


10/28 Administrative Report 

Report on AGARD Meeting 
Resin-Fiber Interface Research 


R . Loewy 
R. Loewy 
J. Diefendorf 


11/4 Administrative Report 

Composite Fatigue Research 

Report on Mechanics of Composites Re 
view Meeting 


R . Loewy 
E . Krempl 

S. Sternstein 


11/11 Administrative Report 

Generic Structural Components 
Computer Aided Design and Analysis 


R* Loewy 
D. Goetschel 
M- Shephard 
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TABLE 

Date 

1983 

11/18 


12/2 

12/9 

1984 

1/20 

1/27 

2/3 

2/10 

2/17 


VI-5 continued 


Topic 


Administrative Report 
Fabrication Technology Experiments 


Resin-Fiber Interface Research 

Administrative Report 

Composites Fatigue Research 

Matrix Characterization and Environ- 
mental Effects 

Administrative Report 
Generic Structural Components 
Computer Aided Design and Analysis 


Administrative Report 

Computer Aided Design and Analysis 

Free-Edge Failures of Laminates 

Report on Gordon Conference 
Fabrication Technology Experiments 


Coupled Bending Torsion in Composite 
Plates 

Composite Fatigue Research 

Unbalanced Laminates 
Quantifying St.-Venant's Principle 

Administrative Report 
Resin-Fiber Interface Research 
Computer Aided Design and Analysis 


Responsible 

Faculty 

R. Loewy 

F. Bundy 
H. Hager up 
V. Paedelt 

J. Diefendorf 

R- Loewy 

E . Krempl 

S, Sternstein 

R. Loewy 

D. Goetschel 
M, Shephard 

R. Loewy 

M. Shephard 
T-L. Sham 

S. Sternstein 

F. Bundy 
H. Hagerup 
V. Paedelt 

I-H. Yang 

E , Krempl 

0. Bauch au 
D. Goetschel 

R. Loewy 
J. Diefendorf 
M. Shephard 
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TABLE VI-5 continued 

Date Topic 

1984 

2/24 Fabrication Technology Experiments 

Composite Fatigue Research 

3/2 Free-Edge Failures of Laminates 
Matrix Characterization 

3/9 Administrative Report 
Unbalanced Laminates 
Quantifying St . -Venant * s Principle 

3/23 Resin-Fiber Interface Research 

Computer Aided Design and Research 

3/30 Free-Edge Failures of Laminates 

Matrix Characterization and Environ- 
mental Effects 

4/6 Fabrication Technology Experiments. 

Composite Fatigue Research 

4/13 Administrative Report 

Report on AGARD Meeting 
Unbalanced Laminates 

4/20 Quantifying St. -Venant *s Principle 
Resin-Fiber Interface Research 

4/27 Computer Aided Design and Analysis 
Free-Edge Failure of Laminates 


Responsible 

Faculty 

F. Btmdy 
H. Hagerup 
V. Paedelt 

E • Krempl 

T-L, Sham 
S. Sternstein 

R. Loewy 
0. Bauchau 

D, Goetschel 

J, Diefendorf 
M. Shephard 

T-L. Sham 

S. Sternstein 

F . Bundy 
H. Hagerup 
V. Paedelt 

E. Krempl 

R. Loewy 
R. Loewy 
O. Bauchau 

D. Goetschel 
J. Diefendorf 

M. Shephard 
T-L. Sham 
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TABLE VI ~ 6 

NASA/AFOSR COMPOSITE MATERIALS AND STRUCTURES PROGRAM ON-SITE VISIT 

December 19 and 20, 1983 

Rensselaer Polytechnic Institute: AGENDA 


Monday, December 19 — JEC 3117 


8:30 

9:15 

10:00 

10:45 

11:00 

11:45 

12:00 

1:00 

1:25 

1:50 

2:15 

3:00 

3:15 

3:30 


4:45 

5:00 

6:30 


Welcome and Program Overview 
Fibers and Interfaces 
Matrix Dominated Properties 
Break 

Characterization of Thermosets 
Discussion 

Lunch (Visitors and Presenters) 

Fatigue in Composites 
Axial -Torsion in Tubes 
Uniaxial Tension at Negative R-Ratios 
Tension-Shear in Tapered Straps 
End-Effect Analysis 
Break 

Fabrication Processes 
Lab Tour 

Interactive Computer Graphics 
Mechanics of Materials 
Materials Characterization 
Thermal Analysis 
Fabrication and Testing 
Discussion 
Adjourn 

Dinner (visitors and Presenters) 


R. Ixjewy 

R. Diefendorf 

S, Sternstein 

B. Wunderlich 


E . Krempl 

* 

D, Elzey 

R. Loewy 

D. Goetschel 

F. Bundy 

M, Shephard 

E. Krempl 

S. Sternstein 
B. Wunderlich 
V. Paedelt 


* 

Graduate Student 
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TABLE 

VI-6 continued 


Tuesday 

, December 20 


8:30 

Finite Element Studies of Fracture 

M, Shephard 

9:15 

Delamination Failures at Edges 

S , Sham 

10:00 

Break 


10:15 

Mechanical Connections in Composite Shafts 

•* 

S, Anderson 

11:00 

Warping Effects in Unbalanced Laminates 

O, Bauchau 

11:45 

Future Plans 

R. Loewy 

12:00 

Lunch (Visitors and Presenters) 


1:00 

Simultaneous Meetings 

NASA/AFOSR (Executive Session) 
Industrial Technical Adbisory Com./RPI 


2:00 

NASA/AFOSR Program Guidance 

Dr. Vosteen 
Dr. Amos 

3:00 

Adjourn 





* 

Graduate Student 
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TABLE VI ~ 7 

NASA/AFOSR COMPOSITE MATERIALS AND STRUCTURES PROGRAM ON-SITE VISIT 

December 19 and 20, 1983 


ATTENDEES 


Dr. Anthony K. Amos AFOSR 

Dr. Christos C. Chamis NASA - Lewis 

lir. Samuel Das tin Grumman Aerospace 

Mr. John T. Quinlivan Boeing 

Mr. Karl V. Stenberg McDonnell Aircraft Co. 

Dr. Louis F. Vo Steen NASA - Headquarters 

Jlr. Max E. Waddoups General Dynamics 
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PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy , Robert G • , Ph.D. 
Wiberley, Stephen E., Ph.D. 

Senior Investigators 

Bauchau, 0. , Ph.D, 
(Structural dynamics, ad- 
vanced composites) *. 

Bundy, F. P. , Ph.D. 

(Physical chemistry and 
structures testing) * 

Diefendorf^, R. J. , Ph.D. 
(Fabrication, resin matrix, 
fiber behavior, interfaces)* 

Feeser^, L. J. , Ph.D. 
(Computer applications and 
graphics , computer-aided 
design, optimization)* 

Goetschel, D. B. , Ph.D. 
(Structural analysis design 
and testing) * 

Hagerup, H. J. , Ph.D. 
(Aerodynamics, configura- 
tion, pilot accomodation, 
flight testing) * 

Krerapl, E. , Dr.Ing. 

(Fatigue studies, failure 
criteria) * 

Sham, T-L. , Ph.D, 

(Fracture mechanics, com- 
posites) * 


•9f 

Fields of Speciality 

^Member of Budget Committee 
tigators 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 


Assistant Professor of Aero 
nautical Engineering 

Research Professor of Mate- 
rials Engineering 


Professor of Materials Engi 
neering 


Professor of Civil Engineering 
Associate Dean, School of Engi 
neering 


Assistant Professor of Me- 
chanical Engineering 


Associate Professor of Aero- 
neering Engineering 


Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Assistant Professor of Me- 
chanical Engineering 


together with Co-Principal Inves- 
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Shephard, M, S, , Ph.D, 
(Computer graphics, finite 
element methods) * 


Associate Director, Center for 
Interactive Computer Graphics 
and Assistant Professor of 
Civil Engineering 


Sternstein^, S, S., Ph.D. William Weightman Walker 

(Failure analysis, matrix Professor of Polymer Engineer' 

behavior, moisture effects)* ing 


V7underlich, B. , Ph.D. Professor of Chemistry 

(Processing science , con- 
stituent material charac- 
teristics) * 


Research Staff 


Manager & Master Technician, Composites Laboratory 
Paedelt, Volker 


Research Associates 

Grebowicz, Janusz, Ph.D Suzuki, Hidematsu, Ph.D. 

Menczel, Joseph, Ph.D. 

Graduate Assistants 


An Deuk, M.S. 


Singh, Sachchida, M.S. 

Burd, Gary, B.S. 


Srinivasan, Krishna, B.Tech. 

Chen , Kuong- j ung , M 

.S. 

Uzoh, Cyprian, B.S. 

Coffenberry, Brian, 

B.S. 

Weidner , Theodore , B.S. 

Gjaja, Niko, M.S. 


Yang, Chung-Ming (Philip) , M.S 

Hu, Tsay-hsin, M.E. 


Yehia, Nabil, M.S, 

Judovits, Lawrence, 

M.S. 

Yerry, Mark, M.E. 

Liu, Shiann-hsing, 

M.S. 

Yurgartis, Steven, M.S. 


Shick, Douglas, M.S. 


Undergraduate Assistants 
Caligiuri, Joseph 
Cas abe 11a, Susan 
Halicki, Paul 
lorizzo, Carol 
Keavney , Thomas 


Seniors 

Lewis, Scott 
Lopez , Matthew 
Schmitt, William 
Williams , Thomas 


Fields of Specialty 

^Member of Budget Committee together with Co-Principal In- 
vestigators 


Undergraduate Assistants 
Basel, Roger 
Botwinik , Sharon 
Cimino, Paul 
De La Fuente, Horacio 
Drake, James 
Father, Richard 
Henderson , Robert 
Hansen, Michael 

Undergraduate Assistants 
Desmarais , Michael 
Egbert, Mark 


- Juniors 

Karkow, Jon 
Kirker, Phillip 
O’Connell, James 
Rajewski, David 
Sawyer, Christina 
Sohn, Kyu Yeon 
Tse, Suet Ting 
Witt, Jay 

- Sophmores 

Kashynski , Stephen 
Morris , Scott 
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